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The effect of atmospheric drag on the motion of an artificial satellite in the gravitational field of an 
oblate earth is evaluated. A spherical exponential model atmosphere is assumed together with the velocity- 
square law in drag acceleration. Atmospheric rotation is neglected. The method of canonical transformation 
is shown to be powerful even in the problem of nonconservative systems like the drag problem. The principal 
effect of the coupling between the nonsphericity of the earth and air drag is shown to be mixed secular terms 
and long-period terms produced by short-period terms of the drag-free problem. The result is complete up 
to the second order in case of secular and mixed secular terms and the first order in case of periodic terms 
if the drag coefficient A* and the coefficient of the second harmonic of the earth potential kz: are considered 
small quantities of the first order. The only exception is with long-period terms in the mean anomaly, for 
which the result is complete to the order A*/k2. The result is complete to the fourth power of the eccentricity 
in all cases, and is valid for all inclination except the critical inclination case. 


INTRODUCTION 


N most actual cases of artificial satellite motion, the 
effect of atmospheric drag is the most important 
next to the effect of oblateness. The usual approach of 


_ combining a theory of the oblateness perturbations, free 


from drag, with a theory of the drag effect on elliptic 
motion is clearly only a makeshift solution of the 
problem. The present paper is concerned with the 
problem of attending to oblateness perturbations and 
the drag effect in a single solution. The approach is 
essentially the familiar one of the method of the vari- 
ation of arbitrary constants. The use of canonical 
variables has particular advantages in the treatment 
of this problem because of the ease with which the 
necessary transformation of variables can be made. 
The development of the theory consists of five stages: 


(i) Derivation of the equations of motion in the 
form expressed with Delaunay variables. 

(ii) Derivation of the equations of motion in the form 
expressed with the canonical variables which are the 
solution of the drag-free problem. 

(iii) Evaluation of drag acceleration in the new form 
of the equations of motion as functions of the doubly 
primed Delaunay variables. 

(iv) Series development of the drag acceleration in 


powers of the eccentricity and in multiples of the mean 
anomaly. 
(v) The integration of the equations of motion. 


EQUATIONS OF MOTION FOR THE DELAUNAY VARIABLES 


The equations of motion in an equatorial Cartesian 
system of coordinates are 


es OU 
ea X 4 j=, 2, oF 
dt? Ox; 

or, in canonical form, 
dé; OF dn; OF 
ain mag ae Sat te oe oe (1) 
dt On; dt 0g; 

if 

£;=2;, j= %), and F=—4> > #7+U. 


The force function U includes the oblateness terms and 
the terms X; represent the drag acceleration. 

It will be assumed that the drag acceleration on the 
instantaneous velocity is directly proportional to the 
cross-section area, the atmospheric density, and the 
square of the velocity V relative to the atmosphere, and 
inversely proportional to the mass. It wili also be 
assumed that the atmospheric density may be repre- 
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sented by a spherical exponential model from the 
perigee height upward. The nonsphericity of the 
atmosphere should be considered in further refinements 
of the theory; this presents no essential difficulties. 

Let w be the angular velocity of the atmosphere, 7 the 
radius vector, and a the osculating semimajor axis of 
the orbit. The expression for X; is then 


X;=—A*C(4;+6,;)V exp(—ar) 


with A* and a constants; C is also a constant, being so 
chosen that C exp(—ar) is of the order unity, 


2 i 
Ve 4|—- — 2 (a%o— x41) + w? (w+), 
r a 
and 
j= + are, 


In the following, the rotation velocity of the atmosphere 
will be neglected. This may be less significant than the 
omission of nonsphericity of the atmosphere. Therefore 
the expressions for X; and V are reduced to 


X;=—AV exp(—ar)é;, A=A*C, 


63=0. 


62= — aX), 


(2) 


Applying to Eqs. (1) the canonical transformation 
(&;,n;) > (L;,1;), we have 


aL; oF dl; oF 2 
ea thy a Oe 
Gh Cgéba| (2 ech Melon ae 2b? 
with 
chia OX}, 
P;=d Mims Qj=h Me (4) 
k OL; k 6 6OL; 


where L;, /; (j=1, 2, 3) are used for the Delaunay vari- 
ables, L, G, H, 1, g, h. The evaluation of P;, QO; leads to 
the results: 


NY 


One 
P;=—AV exp(—ar)p;, pj;=> a 
k 


j 


Onk 
O7= AV exp(—ar)g;, j= 2! fe 
k OL; 
With the use of relations in Kepler motion, it is found 
that 


fi=LL(2a/r)—1], i= 2e sinu+ (2/e)(L2/L1) sinf, 
p2=L», Ope (2/e) sinf, 

ps=Ls, qa=9, 

in which the eccentric anomaly wu, the true anomaly f, 
as well as a, e, r in the usual meaning are functions of 
the Delaunay variables. The relations between P;, Q;, 


and X ; will be proved later as a special case of a general 
canonical transformation. 


(5) 


ATED Gr EEO RS 


EQUATIONS OF MOTION FOR THE CANONICAL VARIABLES} 
WHICH ARE THE SOLUTION OF THE 
DRAG-FREE PROBLEM 


Applying to Eqs. (3) the second canonical trans- 


formation (L;,1;) > (L;",;), we have _ He 
OT py, 2 
as aie Fea Ci 5 Deas A 
dt 4 dt ols" ‘ ‘ (6)| 
Pe = PE a 1 
with a 
OL | 
=) pi Bios Pr 1) 
bo. Ob Se i. 
(7) | 
C= 26 a Q. a | 
z “4 aL; ee e aL; we 


The canonical variables L,’’, 1,’ 


drag-free problem and 
Re* (fe Le l; (Ly Le! Le 

Li Ly La Le lh ole ake 

are supplied by the drag-free theory (Brouwer 1959), | 
provided the inclination is not too near the critical } 
inclination. \ 
Let p;”, gq; be related to P;’’, Q,;’ in the same | 
manner as p;, g; are related to P;, Q; so that 


are the solution of the ) 


OLy 
i=2 ppb él; aA ey Gis gn él; Ww 
P Oly, Oly 
qi =2 Di aL, ee Vises aL, Fe 


Then p;”, q;’” may be written 


bi =p;+5p;, 9; =qit6q;, 


with 

0(1,— 1") Oly 
7 86i= 2 Pi ee Cian pt 
k él,” Olen ( ) { 
fae it 
Al A(Lx— Ly") i 

éqj= Pri ae dk 5 
k OL; k aL ;” 


6p;, 6g; being of the order of the oblateness perturba- 
tions which are considered to be small quantities of the 
first order. Assuming the drag coefficient also a first- 
order small quantity at the most, Pir Qi present in Eqs. 
(8) may be replaced by the expressions (5) with change 
of (L;,l;) to (L,’",0;’") if effects up to the second order 
are considered, as is the case. 

A general proof of canonical transformations of the 
type Eqs. (4), (7) is easily given. Let 


dx; dy; oF 
= h —=—-——V;, 


oe 5 NWN 
dt Vi dt Ox; 


j=l, 


be given, in which X;, Y; are arbitrary functions of the 


| 
| 
| 
| 
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‘) variables «;, y;, and t. Let 


Urange 


4 
) 


4 


f the 


ame 


i 
| 
; 
i 
: 


=H ;(Le Ye), Ve=Vi(Xe' Ve’) 
| be a canonical transformation and 


dx; oF* dy; 0F* 


mz jo j 


Oy; dt Ox 


F (23,93) = F*(x;',y;') 


dt 


with 
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In the same way, 


yee 
oe) 
Xp 
EVALUATION OF 69;, 6q; 
If (p;)", (q;)” stand for p;, g; defined by Eqs. (5) 
with the elements L;, /; replaced by L,’’, 1,;’’, Eqs. (8) 
may be written 


: : O(k—lk””) 
be the transformed equations. nee pS yn 
With use of the original equations we have y 7a al,” a ed 
(8’) 
dee; Oy; 5 Ox; OF* OY; Ox; Ol; 0( p— Ly”) 
j= tz (x +7). 89)= EC)" 
dt OYy! a Or 7 OF, 1 Oyy! dy, k oe OL,” 
Tn the left-hand member, the substitution The evaluation of 6;, 6g; with the use of the above 
noe Oe ie ee ih! equations and the solution for L;, /; in the drag-free 
Al ie aca ) theory is straightforward but laborious. By taking 
di & \Ox' dt dy! dt advantage of properties of canonical transformations 
f i the work can be greatly reduced. 
dy =¥/( dy; dete mn O95 Wt ) The expressions for (p;)’’, (qj) are, in the original 
di \dx! di dy,' di form, < é 
yields Mt Oe Om 
day! dy;! ax, (p;) =: & ) COE = & pel 
E (Lave F——+in'o9'— = vals" aly 
: ‘a gi e Introducing the determining functions S$, S;* which 
- since all the Lagrange brackets are zero except were defined in the drag-free theory, we have 
[¥p',¥p' J=+1. Rie *) 
_ We have then oe ae ! 
dxy) oF* * yi) 3 (9) 
Fhe Ras Tay) ? — Ae ———e 
ay, ae aye! ie Gh serra 
) 
Be nally aye where it is permitted to replace J; by 7,’ in S; and 1,’ by 
eet -E (x= : ). 1,’ in S;*. Therefore, if partial derivatives with respect 
Ay," to 1,’’ or L;’’ are expressed by subscripts, 
Vp p V Pp 
Ons ’f 9(S11;"+S11;’*) One’P = 0(S11;7+S11;’*) 
ip EE | - pepe || 
k i Ol,” a oe a aaa Bg dl; 
BS cota pe dE! a(Suy" +S") a 
k z ol; D 0&,” OL,” ony” OL,” 
Oni” O(S11;”+S11;’*) OE)” A (S11;77+S11;’"*) Any” 
1h % Ss g" : J - J ~ J F] Pp 
kG OL, »? dé, Ol; f Ony” 01," 
(S117 +S11;"*) Oni” 0&,"" Oni” OE,” 
=DDEEH" |- | 
k ip 0&,” ole’ OL," OL,’ EA 
0(S11;7+S11;'")F One’ On” ~— ny"! Any” 
+D 0D &’——__ | - =| 
k i p Ony” Ol” OL,” OL,’ Ol," 
, O(S1yi"FS11;"*) oe 2 0 (S143 +S11;'"*) . 
= ee {71 F s+ ~ et stn 
‘ae ae ony” 
= pears tees”) : 
i Oe 
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since all the Poisson brackets are zero except 
{ni ,é2} es ils 
In an analogous way, we have 


0(Sinj”+S11;"*) 
Wee ee, 
i 0&; 


Introduction of the operator D defined by 


(10) 


Sete - 
i 0&; 


leads to simple expressions for 6p; and 6q;: 


re OS; - OSy aka 1) 
bpi= ( ~ ). 6q;= ( ). 1 
4 qi oL Olle. 


ol’ 
With use of a property of canonical transformations, it 
is easily verified that 


0Si* 
ol,” 


We thus find 
WwW» 


AGN DR GF EVO IRD 


OL,” Onn’ 
Die -> &,/ = —-> &,/" aes (p3)”, 
k Eq _k Olea | 
and (12) i 
él; yr Onn” 5 i} 
DN fa Se = +h a =+(q;)”. 
k 0&,”” Cy 


An immediate consequence is 
DEL) = 0 , 


or, if J stands for the orbital inclination, 


DI" =0, 
By the definition of D, we also have 
Wr, 
or, 
Dr —=De"= 


or (q;)". 


5p1=3 ; 


Spo= 


GAaND 


2Ro 


ee 1662+ 156°) (1— 56”) 


wR 2 a Si BS) 
= | (-14+30) (S441) sinas-+-(“—“e) | 
2L'G r Yr DBS i) 


wh @ 
eL’G 


fi al 
Te = ) sin(ae+4) I+ 
ha FILS 
OMS OF 
jeslena 
US Dy r 
wks 1 
ad 
eL'G 
qa 
=| (-1+30(=n 
vba 


es ate -) sine) |= 


2k 


58 
a. “+1) sina (———o 
2 


6q2= <= 
BACHE? 


ae 2B 


Bs? 


eL?G? ¢ 


a 


o 1+36) (<a 


| 


=| (- 14309 (— pee “+4) sin f 


‘lk a a 5 a a’ a 
()|-" 451-2) | H(; ) i) een) ; 


WEIS 1 1 
( —#)|- cos(2g-+2f)+e bose nT a cos(2e-+3,)| 


1 1 
AR. cos2g-++-— cos(2et 
8 a 
Ga 
ee, sin2g 
Te r 


2 


PS “+4) sinf 


a a 
sin(2e¢+ f)+ (Zr+2+2) sin(ag+3)) | 
ie ip 


: A = 1167) = 10081 -s#)>]| ( i \e sin2g+sin(2g+ f)— sing) 
r 


e a 


eet sin2g 
r Yr 


(14) 


+( 0 )\[(—=» ~“+2) sinh +( “r+"42) sin(2e+3)]} 


oY r 


2Ro 


ae 3 (o ae 1) sing (-—“#)| ——“1) sin(2¢+ f) 


(13a) | 


(13b) | 


8 being the latitude of the satellite. The meaning of | 
double primes in Eqs. (13) is the same as that of (p;)’” |) 
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a ar 1 pk 1 
+(Sr+-+-) sin(2e-+39) |} "(34150 207) +e(6inf—sinu)+—(0—w sinf| 


wR 1 1 1 wkoyp 1 11 
Tes g 150| - sin(2g+2/) a Saab a sinQe+3/) |] +e +2064(1—56? | sin2g 


2k 


“—- (1—336) — 2564(1—56?)-!— 506° (1— se) Dsin(2g+-f)-+sin( 2g—f)] 


who 1 37. 55 
<= | <4 — 330 nt (1-—# ) + (—2597-+35)0"(1— 58) SOG or —se)-] 
eG L8 See 3 


X[sin(2g+ f)—sin(2g— f) ], 


2R 


1 1 1 1 
— -6q3= o—a2 (f—))+e(sin f—sinw)+—(1—7) hoe eine Dee ete dr sin(ae+39)] 
e 


2k P11 
— f+ 108 —56°)1-++2564(1— soy] sin(2g+f), 


| where double primes are omitted for simplicity of printing, and 


0=H"/G" n=G"/L". 


The contributions of the fourth and third harmonics are obtained with use of A,S;* and A3S;* in the drag-free 


theory. The result is, double primes being omitted, 


| 5 pw ka 
Adpm— tle 36?— 864(1— 56?) ]Le cos2g+4 cos(2g+f) ], 


1 wp Az 


2 
A;:6p2=- — igen snl sing-+- sin(e+ |; 
4 G ko é 


5 pe Rs a , P 
A.6qi1= —- [1-39 801s)" e( 1-3") sindg-sin(2g-+/)—sin(2e—)} 
. 6 eL’G ke r 
8 snl e(1- =) (e—f) 
ANG = cosg—cos(g— | 
ce ae, 2 
we? Ra we Rg 5 50 100 
soe — Par = ~o(1— 5@°)- 2] snag =e —— 190) +91 54+ — 40150) 
G* Rol3 G+ ke 1 3 3 


(14’) 
we kg 70 
X [sin (2g+ f)+sin(2g— ine lay ae (i= se)+(~ <i —— Jo —56) 


100 
+= (7 Dom(1— 50") [sin e+) sine NI 


1 pw Az ofsinl 
dap 


one =|? cose cos(g— nj conte tf); 


é 


2 


Lu 4 16 40 
As5q3=5— etl 1+ P5140 — 50] sin(2g+/f), 
. Gtke 3 3 


1 yo A3.9 6 
A36q3a=- — ——- —— cos(g+f). 
2G ky sinl 
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| 
It should be noted that short-period terms, when SERIES DEVELOPMENT IN POWERS OF THE ECCENTRICITY, 
operated on by D, produce secular and long-period AND IN MULTIPLES OF THE MEAN ANOMALY 


terms. A typical example is the derivation of 61: The equations of motion, Eqs. (6), may Delmieen ink 


| ist 
OS; pho toeS q’? the form yy 
cat 2 ae a. 
hoa —|( alee )( n an, dL ,!'/dt=—AV exp(—ar)p; 
—AV" exp(—ar’")5p; | 
Ones He { xp( )6p;, (15) 
+(---#)— cos(2¢"+2/")] dl,!'/dt= — (0F**/0L;’)+AV exp(—ar)q; eA | 
Dn Died fea 


+AV” exp(—ar”)éq;. 
oi only short-period terms, while 51, which is derived Owing to the presence of the factor V exp(—ar) or| 
Mf 


Fy V" exp(—ar”) in the right-hand members of the} ; 
Spem Diebabalc) equations, a series development in terms of the eccen-|\j 
has both constant and long-period terms in addition to tricity and the mean anomaly is inevitable before] 


short-period terms, that is, integration. The terms i 
Whey 1 3 \fa"% V exp(—ar)p;, V exp(—ar)q; | 
[Ps lconse= =15; sf ic + )| i 4 
i Died r const 


may be developed in terms of e, / at first, and finally in |j 

7k te terms of e’”’, /’’ with use of the solutions of the drag-free |) 

[1 liong-period = ey si ‘(-- problem, while the terms i 

ng-period = a) 
LEONG Ge) 


V"” exp(—ar’)ip;, V” exp(—ar’’)6q; 
"4 
«|= eos? ‘| -cos2g, are ready to be expanded in terms of Ae 
7/4 


const If now : : | 
x=(r/a)—1, «’=(r'"/a")—-1 (16) } 

This fact is one important result of the coupling effect 4 

between drag and oblateness. are introduced, 


V=(u/a)'[2(a/r)—1}) 
= (u/a)*[(1—x)/(1+«) }! 
= (u/a)* (1-430? — 303 + gat — 34°), 
exp(—ar) =exp[—aa(1+<) ] 
=exp(—aa) exp(—aaz) 


aa aia atot aba \ . 
exp(—aa) (1 aarti at, (17) | 
2 24 120 y 


1 1 ike! | 1 1 
V exp(—ar)= (n/a) exp(—aa)|1— (Ian) (<a raat (<4 oat atatt oat) 


feels 1 ft if Saks 1 1 1 1 : 
ae (Gant etat cat +— tat at — ( —+-aa+—o?a?-+—o'%a+—otat+ oe) (18) — 

come ee 4 6 24 Sus “ 12 24 1205 5 
By the definition, x is a variable of the order e. Convergence of the above expansions depends on the magnitude 
of aae. Tentatively the expansion is carried to the fifth power of aax in order to make the solution complete to 
the fourth order. If age is large enough to endanger the convergence, a remedy may be to use 


2 — 


ane 
6 24 120 
in place of the Taylor expansion used above. The necessary change is: 


A 20a? A 30°03 A sota* A saa” 
exp(—ar)=exp(—aa)| Ay Aart a] 


(a7) | 


— 


at — A jai, j=9, i 3 5; 


and may be done at the very end of the development of the theory since powers of a are preserved without any — 
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ny change throughout the theory. The numerical values of A; may be best determined when a value of aae is given 


| in an actual case. 
Although V exp(—ar)p;, V exp(—ar)q; should be considered as functions of Delaunay variables, it is con- 


. . venient to express them as functions of a, x, e, f, and J. The explicit expressions are 


. V exp(—ar)pi=V exp(—ar)L[2(a/r)—1] 


9 1 11 9 3 1 
Re | =p exp(—a| 1— (3-++aa)a+ (<-+Sa0+ ota )at— ( +—aa-+-0a?—+- a) 8 
2 2 Die NPs v 6 


. 

1 1 eYoNoet| 11 
Se sat ) -(—+- aa “ato “eels wou 155 ote! at] 
j 


| 
| V exp(—ar)p2=V exp(—ar)G 


Ht exp(—aoh| t= (Ataaye+(~ a - (5+ +-aa-+ “ta <a)? 


| a | 1 1 Pale 1 ieee (aes 
| ~§ + aa-+-ot-+-atu-+— ata) oa (+ -00-+-atat+—ate+— ata +—ata') | 
A aera 24 Sie. ene 12 24 120 


V exp(—ar)p;=V exp(—ar)H 
| =[V exp(—ar)p.10, 


1} 2 
_ V exp(—ar)qi=V exp(—ar) (2 sinu+— sn) (19) 
€ 


b\? 1 il el 1 1 
= (*) exp(—aa)| 1 (1+aa)e+ (<ta0+-o'at ) (+00 tt-baatat) 8 
| a 2 2 DOM ae ING 


aa | 1 1 Tmo 1 1 1 
| + (<a -oPat + -oPat + — rat )at—( +-aa+-—oa?a?+— 08 a? +— oat ata) 
OMe 4 6 24 8 8 4 119) 24 120 

| Ten es 1 
| x[-(cet#) sinf-+2e( 14-2) sin 
1\) e€ 2 8 2 

2 
V exp(—ar)q2=V exp(—ar) (-- sn) 

e 


uy? 2 1 1 ft 1 1 
=— (*) exp (—aa)— snj|1 — (i+-aa)x+ (<-taa+-otat)at— (<tcatnatat nate! Jat 
a @ 2 2 22 yy 6 


Ee i 1 1 he 1 1 1 1 
a (<+-a0-}-atot-+-otal+— ota ) ee (<-+-aat-ot0+ ofa 4-——of att ota | 
moe, 4 6 24 8.980 oA 12 24 120 


V exp(—ar)q;=0. 


If ¥ stands for any of these functions, we have x= (r/a)—1=—e cosu, 


| ! ayy” ayy" ayy" dx/de=—cosf, dx/dl= (e/n) sinf, 
vay'+ (—) ba (~) bx-- (~) be therefore, 
. da Ox de 
ap+(— ). 62. (20) while, 


+(=) id 
. 8f= (Af/de)be+ (Af/al)al 
The variations 6a, de, 62 are directly available in the a 2 2 Y inf-de-t a CA (€2/2) 
result of the drag-free theory, and the variations 6x, =1*(2-+e cosf) sinf-de+n“L1+ ( 
6f should be expressed in terms.of 6e, 61. +2¢ cosf+ (e/2) cos2f]él. 


dx= (dx/de)de+ (dx/ 01) dl 
= —cosfde+ (e/n) sinfél, 
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By the drag-free theory, on the other hand, we have 


2k es: a Seo Ne 
so-—| (- +4) ~—9+)+(-—#)— cos(2e +21) 
a aed a DI TS A a 


be=b1ge+5sne, 


ky 1 Az sinI | 
bye —[1—116°—4064(1—56?)- "-- “11 3@°— 864(1— 562) ] | el So 5 —— —— sing, 
a n? 8 4 ko a 
ko 143 a She a 
ine=— a lee (=-r*) a (---#) (=-r+) cos(2g+2/f) 
we ie? r hay) rs 
1 


33 
ae 0 )3e cos(2g+ f)+e cos(2g+3f) ]¢, 
SiN) 2 


61=bigl +6sn1, 


61,f= = 519, 


7? tanl 


i 
inl = 4 sinI[3 cos(2g-+2)+-3e cos(2e+f) +e cos(2g+3f)], 
an 


51=8,1+8i9l+6en/, 


3 
v= +0) 
(Pips. DY 
ko yf 1 Azo sind 
sua —[1—116°— 4064(1—562)- te woe 36°— 864(1— se)-4}} sin2g——- —-——y cosg, 
anl8 12 ke? ko | ae 


@en Dee 


ke Sere e a 
6 eo= -=| (-—+ #) (Ses +1) sin f-+-— (--)[ ot) sin(2g+f) 
r r ey 


GC Gin ty 
+ (Gr+-+-) sin eetsp)| ; 
r TAS 
With these formulas, an explicit expression for 6x is 


6x= —cosf(Sige+5sn€) + (e/n) sin f (5.1 +619) +6,;/). 
The following is the result: 


6x=d4+6nx+dmr4, 


kp 
vos sin f 6,J= 3— —-+0)< sinf-nt, 
n ant 


é 
dux=—cosf d1,¢+- sin fbi) 
UT] 


=— “{1- 116°— 4064(1—56?)—1 ee 2 p)—1 ae 
a (1-509) [13 —801(1-5)-"T fe cos(2¢ +f)" — sin(e tf) 
2 
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ae (21’) 
durx= —cosf 6 ro, sin f dsal 


3 


SS0LC 


} ers \ (1 fa lia da 
i HO) {-(Sat—1 )fcos(2e-+ f)-+eost2e+3/1]-5(<a+"-+1) cos2g 
Dine 2 e\r 2\7r r 


=n) age (= Zs “+1 )(1—cos2/)] 


Pa? 


9 


a 7 5 ete 1/@ Le 5 
+(> +-—-) cos(2g-+ N- (= ++") cos(2e+4f)|). 


r r r 
Since f is to be considered as a function of e and /, we have in Eq. (20) 
(oy/af)"5f=Loy (f)/ol ]"sl+ Loy (f)/de]'Se, (22) 


the latter part of which is absorbed in (d/de)’’5e. It is understood that all the variables in the expressions for 


da, de, 61, 6x are doubly primed. 
The variations, 6a, 6e, 61, 6x should be also developed in powers of «. The final resale for these expressions is 
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+[-S5 set (s-e) = (28 e420 sx cosag+| =e (s— ats (4—3e?)x? 
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It will be noticed that 5,,/ appears in the calculation of both / and g; when / and g are combined to /+-g, th 
result is complete to the order e!, though the above expression for 5,,/ has a divisor e. 
With use of Eqs. (20), (22), we may write the result of expansion in the following way 
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where, for instance, (p;) means [V exp(—ar)p, |’ and (p;)2 means 
{(9/da)LV exp(—ar)p; da}. 


The meaning of the other symbols will be clear without confusion. After a lengthy calculation, we find: 
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In the foregoing equations, all the variables in the right-hand members should be understood to be doubly primed. 
The development is complete to e*. The divisor e” in the right-hand members of the equations for dl’’/dt and dg’’/dt 
will be canceled in d(/’’+-g’’)/dt; the divisor sin/J in the right-hand members of the equations for dg’’/dt and dh’’/dt 
will be also canceled in d(g’”+h’’)/dt. As to the divisor, 1—56?, the situation is different and the expansion may 
become invalid when the inclination is near the critical inclination, where the solution of the critical inclination 

(Hori 1960) may be used as a drag-free solution without any change of the method presented in this paper. 


INTEGRATION 


The‘method of successive approximation is applied. The zeroth approximation (A=0) corresponds to the 
solution of the drag-free problem: 
Lj'=const, 1,’=—(0F**/dL;')\t+ const. 
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the equations for dL ,”/dt are integrated, then the equations for d/;'"/dt. In the terms — 0F**/dL,'’, since they i) f 
not contain the drag-coefficient A, the first-approximation results of L,’’ should’ be substituted except the par} 
that contain k,”. i 

With use of the equations of dL,/’/dt with a, e, I constants and J, g linear functions of the time in the right-han} 
members of them, the first approximation solution for L,/’ may be obtained after integration. By the drag-fre}} 
theory, we have 
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where Ly”, Go”, and Ho” are three integration constants. Before proceeding to integration of 1’, g’”, and h’’, the 
term —0F**/dL", —dF**/dG", —dF**/dH” should be evaluated with substitution of the above results in L’’, 
G”, H” included in them through a, e, I, n. 

By the drag-free theory, we have 
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pas the variables in the right-hand members should be understood as doubly primed, and functions of L”, 
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In thefexpression for 7, the terms of the order A*k» should be included as well as the terms of the order A* in AL, 


while in the expressions for the others, only the terms of the order A* are needed since the factor k2 is multiplied. 
With the results for AL, AG, AH, we obtain 
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After integration we have as the first approximation solution for 1’, g’”’, h’’, 
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1y’’ being a constant of integration. 


As can be seen by the expression for nf’ (—3AL)dt, the accuracy of the long-period terms in /” is limited to the 
order A*/k:. This is due to the presence of the term m without the drag coefficient A nor the coefficient of the 
oblateness 2 in the right-hand member of the equation for d/’’/dt through the term — 0F**/0L”’. In order to obtain 
long-period terms to the order A%*, ks” terms are needed for long-period terms of L in the drag-free theory, which 
are not available for the moment. In the drag-free problem, the situation may be the same if integration of long- 
period terms is performed with respect to ¢/ as is the case in the present paper, but the difficulty was avoided by 
performing the integration with respect to g with the help of the canonical transformation method (Brouwer 
1959; Garfinkel 1959). This does not apply to the drag problem owing to the simple fact that the system is not 
conservative, although a method of canonical transformation can simplify the intermediate calculation needed. 
In the cases of g’” and /’, long-period terms are complete up to the order A* since the terms —dF**/dG” 
and —0F**/0H” have the oblateness coefficient k2. The following is the result of the integration. 
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These are the result of the first approximation. The 
second and higher approximations will give terms 
factored by A**/k2, A**/k.?, --- for the variables L’”’, 
G”, H”, g", k” and terms factored by A**/k2, A**/k3, 

- for long-period terms in /’’. The convergence of the 
solution then depends on the magnitude of the ratio 
A*/ks, which is smaller than unity if effects of the drag 
can be assumed smaller than those of the oblateness. 

The expressions for L”, G’, H”, I’, g’’, h’’ obtained 
should be substituted for the corresponding symbols in 


the drag-free solution (Brouwer 1959) to give the 
solution with drag. The six constants of integration of 
this solution are Lp”, Go”, Ho’, lo’, go’, ho’. 

The completely literal character of the theory makes 
it applicable to any actual case of satellite motion, but 
it is important to be aware of certain limitations that | 
may effect the effectiveness of the theory in particular 
cases. In this connection the assumptions made for the 
earth’s atmosphere are particularly relevant: (1) for 
the atmospheric density we have assumed the ex- 
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pression poexp(—ar), po being a constant, which 
implies spherical symmetry of the atmosphere; (2) the 
rotation of the atmosphere was ignored. 

Even if in specific applications these simplifying 
assumptions appear justified, there remain limitations 
introduced by series developments. For the formula for 
the atmospheric density we introduced into the theory 
a Taylor expansion 


Pat) — 1-255" —En"---, 


in which x takes the value —aae at perigee, +-aae at 
apogee. For cases with large values of aae the con- 
vergence is insufficient. In the section on Series Develop- 
ments we recommended a modification available if 
actual numerical values for a, e are introduced. Use of 
the coefficients Ao, A1:-- will improve the convergence. 

Finally, the effectiveness of the theory depends on 


‘the value of A*/k». For satellites with low perigee 


heights the large value of this ratio is unfavorable to 
the convergence. 
The presence of k» in the denominator in coefficients, 


_ such as in long-period terms in the mean anomaly in 


the first order, raises the question what becomes of 
such terms in the limiting case ky — 0; 1.e., if the oblate- 
ness is neglected. It is readily seen that in this limiting 
case the relevant terms become linear functions of the 


time or constants, in view of the fact that 


sin sin 
(g,2g)=  (go,2go)=const. 
Ss 


cos 


lim 
ko, ks, Az,0 0 co 
: sin sin 
lim if (g,2g)dt=t  (go,2go)+-const. 
ko, ks, Az,0—0 cos cos 


Such terms may therefore be absorbed into the 
constants of integration. 

In spite of its limitations, the theory has a consider- 
able area of applicability. Moreover, various refine- 
ments and adaptations can be introduced later when 
the need for them arises. The manner in which oblate- 
ness and drag are coupled to produce mixed secular 
terms and long-period terms from short-period terms 
in the drag-free theory appears to be of considerable 
theoretical interest. 
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The effect of the most important longitude-dependent term of the geopotential on the motion of artificial 
satellites is expressed by simple formulae. Its numerical coefficient 6 can be interpreted as the ellipticity of 
the earth’s equator; the phase constant Ax gives the geographical longitude of the semimajor axis. Based 
on photoreduced Baker-Nunn observations of the satellites 1959 a1 and 1959 7, the values of these constants 


were determined to be 
B= (3.2140.29) x 10-5 


VOLUME 66, 


and \x= 


NUMBER 5 


—33°15-0°53. 


INTRODUCTION 


A Re primary mission of satellite 1959 al (Vanguard 
II) was to scan the earth’s cloud cover, and that 
of 1959 7 (Vanguard III) to measure the earth’s mag- 
netic field. These experiments, however, could be 
conducted only for a limited period, as the transmitter 
of Vanguard II ceased operating 27 days after launching, 
and that of Vanguard III 85 days after launching. 
Since then the Smithsonian Astrophysical Observatory 
has had the responsibility for tracking both objects, 
and has obtained many thousands of optical obser- 
vations of them by means of the Baker-Nunn camera 
stations and the voluntary Moonwatch teams. Because 
of their high quality, these observations became very 
valuable for different research projects, such as the 
study of the upper atmosphere, the determination of 
the earth’s gravitational field, and other investigations 
now under way. For our present purpose, the deter- 
mination of the coefficient and phase of the first longi- 
tude-dependent spherical harmonic in the geopotential, 
only the most accurate observations can be used. These 
are without question the positions obtained by the 
Photoreduction Section at the Cambridge headquarters 
of the Smithsonian Astrophysical Observatory from the 
Baker-Nunn photographs. For simplicity, we call these 
positions photoreduced observations as opposed'to the 
field-reduced observations; the latter are the less 
accurate positions determined in the field, i.e., at any 
one of the 12 Baker-Nunn camera stations, immediately 
after a pass of the satellite over the station. 


GEOMETRICAL INTERPRETATION OF THE GEOPOTENTIAL 


Let X, Y, Z be the coordinate system of the principal 
axes of inertia with the origin at the center of mass, and 
ASBSC the moments of inertia about these axes. 
Then the gravitational potential of a body can be 
written in the form 


de eee Hh +0 +0(- ) ay 


iF 2m 


where u=fm. is the product of the constant of gravi- 


tation and the mass of the earth, and 


[=r?(AX°+BY?+C2) 


{ 


is the moment of inertia about the line joining the| 
center of mass to the point considered. It should bel 
noted that, according to the convention in theoretical) 
physics, we are dealing here with a force function, which} 


is the negative of potential energy per unit mass.| 


Expression (1), sometimes referred to as MacCullagh’ S|, 
formula (Jeffreys 1959), is usually expanded into a} 
series of spherical harmonics as follows: 


Mt Ge\* 
v="| 1-7.) Pa(sin®) 
r a } 
=e, 1.0(= A P,2) (sin®) costa. (2) 
E 


Here @ denotes the geocentric latitude; A the longitude 

reckoned from the principal axis X; the P’s are 
Legendre functions; and @, stands for the mean equa- | 
torial radius of the earth; furthermore, 


2C— (A+B) 
b Pmen TES eIa is 


2m Ae 


A-—B 
and. Jip@e= : mi 
4m G2 


Because all gravity measurements are made with 
respect to a frame rotating with the earth, the geo- 
potential W is the sum of the gravitational potential U 
and the potential of the centrifugal force; that is, 


V=U+W(+Y), iee) 


where y is the angular velocity of the earth’s rotation. 

To express the constants J; and J;®) by geometrical 
quantities related to the shape of the earth, and by a 
quantity related to y, we assume that the surface of the 
earth is very nearly an equipotential surface of the 
geopotential defined by Eqs. (2) and (3). It is easy to 
see that these equipotential surfaces for r~d, are 
approximately ellipsoids with three unequal axes. Let 
a> b.>ce be the axes of the equipotential ellipsoid that 
fits best the real surface of the earth, and let the earth’s 
oblateness and the ellipticity of its equator be defined 
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by the relations 


Ce= G,(1—Y) and b-=a-(1—8), 
where 
a= 3 (a+b) “4 


(The above definition of the oblateness for a triaxial 
ellipsoid was suggested to the author by William Kaula.) 
Also, we introduce the small dimensionless quantity 


a=~/(u/a.'), (4) 


which is the ratio of the squares of the earth’s angular 
' velocity and the Keplerian mean motion of a fictitious 
satellite with semimajor axis a=d,. It can be shown 
then that the quantities a, 8, y determine the coefficients 
in the gravitational potential as 


De Way oue 
gee ay— sey (5) 
Sie Oma 
and 
Jo =— 48. 6) 


The value of the constant a can be determined with 
high accuracy. A sidereal day equals 86 164.10 sec; 
therefore, Y=7.292115X10-* sec!. Then, adopting 
from O’Keefe, Eckels, and Squires (1959) the value 


u=3.98618 X 105 km? sec? 


and the value 
G.= 6.378388 X 108 km 


of the international ellipsoid, we obtain from relation (4) 
a= 3.46165 X 10-. 


The most accurate value of the coefficient Js is probably 
that of Kozai (1961), who has found that 


J ,= (1.08219+0.00002) x 10~. 


With these data, Eq. (5) gives for the oblateness of 
the earth 
y= 3.35227 X10-*= 1/298.305. 


The data on the ellipticity of the equator are rather 
contradictory in the literature. They were determined 
by analyzing the measures of gravity over extended 
areas of the earth, and seem to be influenced signifi- 
cantly by the distribution of the observational material. 
For instance, from the treatise of Heiskanen and Meinesz 
(1958) we list the data given in Table I, where x is 
the geographic longitude of the direction of the prin- 
cipal axis X, so that for any longitude \ we have 


A=dA—-Ax. 
According to Jeffreys (1959), however, 
B=(2.540.9)X10-5, Ax=0°. 


In his textbook on celestial mechanics, Subbotin (1949) 
adopted the value 


8=1/30000=3.3X10-°. 
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Taste I. 
8 Ax 
Helmert, 1915 3.6X10-5 —i7° 
Heiskanen, 1924 5.4xX107-> +18° 
Heiskanen, 1928 ere SG I te 0° 
Niskanen, 1945 4.6X10- —4° 
Uotila, 1957 77s Ba Ut Fa —6° 


In his detailed study, Jongolovich (1957) derived the 
equivalent of the values 


8=3.57X10-, 


Ax=—7°7. 


From a recent paper by Kaula (1959) one can infer 
the values 


B=3.7X10-§ and Ax=—20°8. 
EFFECT OF THE ELLIPTICITY OF THE EARTH’S 
EQUATOR ON THE MOTION OF SATELLITES 


Several investigators have pointed out that the 
ellipticity of the earth’s equator might have a detectable 
effect on the motion of artificial satellites, and therefore 
a careful analysis of accurate satellite observations 
could lend itself to a better determination of the 
constants 6 and dx. In this connection let us refer to 
the papers by O’Keefe and Batchlor (1957), Robe 
(1959), Cook (1960), Sehnal (1960), and Musen (1960). 

To derive the perturbations caused by the presence 
of the sectorial harmonic in question, we can confine 
ourselves to a very simple theory. The disturbing 
function, which is 


Bu fa? 
R=-— ~) cos*® cos2A, 
ZG t 


must first be transformed into a function of the orbital 
elements a or n, e, J, M, w, and © of the satellite. Let 


v=the true anomaly, 
u=v-+w the argument of latitude, 
w=the projection of # on the equator, that is, the 
longitude of the satellite reckoned from the ascend- 
ing node of its orbit, 
§¢ and 6x=the mean sidereal time at Greenwich and at 
the intersection of the principal axis X with the 
surface of the earth. 


Note that 


and 
Ax=Ox—9¢. 


The angles ®, w are related to the angles J, u by the 
formulae 
cos® cosw=cosu 
and 
cos® sinw=cos/ sinu, 


in accordance with spherical trigonometry. These equa- 
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tions can be written in the complex form 


cos® exp(iw)=cosu-+7 cos/ sinu 
=cos*4I exp(iu)+sin?3/ exp(—iu), 


which immediately yields 


cos® exp(iA)=cos*4I expli(u—0x+2) | 
+sin?4J exp[—i(w+6x—2) ]. 
But 
cos*® cos2A= Re{[cos@ exp(zA) P}, 


so that introducing the abbreviations 


c=cosl, s=sinl, #=p/a', 


we obtain the disturbing function in the form 


R=}6a2n*(a/r)*{ 2s? cos2(6x—Q) 
+ (1+c)? cos2(u—6x+Q)+ (1—c)? cos2(u+6x—)}. 


General formulae for the disturbing function of any 
tesseral harmonic can be found in a paper by Groves 
(1960). 

In our problem it is not necessary to carry out the 
development of this function into a Fourier series in 
terms of the mean anomaly M, because the contribution 
of the short-periodic terms to the perturbations turns 
out to be negligibly small. As the Fourier series of the 
functions 


(a/r)? cos2v. and (a/r)* sin2v 


do not contain constant terms, and the constant term 
in (a/r)* is (1—e?)-3, the “long-periodic” part of R 
becomes simply 


R= (36)a2n?(1—&)-3s? cos2(@x—2). (7) 
The perturbations caused by the oblateness of the earth 
concern us here only to the extent that Q is not constant, 
but rather a linear function of time: 


Q= +01 (t— to), 


6x—Q=6x9—Q+ (6,;—Q)) (t—to). : 


Substituting function (7) into the differential equa- 
tions of the orbital elements [see, for example, Moulton 
(1959) | and integrating with respect to time, we obtain 
at once 


so that 


na? (3— 5c?) 


o= sin2(@x—{), 
8(6:—-1) p* 
NEC 
6Q=6—— sin2(6x—Q), 
4(6:—Q1) P? 
(8) 
naes 
él =8———— c0s2(6x—2), 
4(0:—-Q1) P? 
3na?2(1—e?)3s? 
M = 6———- sin2(6x—Q). 
8 (01-21) p* 
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| 
The semimajor axis a and the eccentricity e do not | 
undergo long-periodic perturbations due to the ellip- ; 
ticity of the equator. The question now is how to use | 
these formulae for the determination of the constants 8 |) 
and x, and whether the accuracy of the observations | 
made with the Baker-Nunn cameras is sufficient for \ 
such a determination. | 


COMPUTATION OF THE CONSTANTS 3 AND Ax 


In the Differential Orbit Improvement Program of | 
the Smithsonian Astrophysical Observatory developed | 
for the IBM 704 computer by G. Veis and Ch. H. } 
Moore, the orbital elements w, &, I, e, M can be repre- - 
sented as polynomials up to ‘the seventh degree, plus | 
sine terms and exponential functions of the time /. As | 
many as 23 of the data figuring in such a representation | 
can be varied in an iterative least-squares solution, but 
we usually limit ourselves to a much smaller number of | 
unknowns in order to obtain meaningful results. It — 
takes some experience to decide properly which quan- | 
tities should be varied in a particular application. 

The numerical material for the present investigation 
was obtained in December, 1960 by several consecutive - 
runs of this program, as applied to 187 photoreduced_ : 
Baker-Nunn observations of satellite 1959 a1 covering 
the interval April 6 through 25, 1960; and to 216 obser- 
vations of satellite 1959 » covering the interval May 8 
through 27, 1960. The use of two satellites is not a 
theoretical necessity, but a precautionary measure in 
order to give two independent determinations of the 
constants in question. While it is not possible here to 
furnish all the details of these orbit computations, their 
important features can be summarized as follows. 
Because the phase at which the trigonometric terms (8) 
are “observed” at a particular station varies only 
slightly from day to day, it is desirable to use a long 
interval of time for the computation of the constants 
6 and \x. However, even in the case of such relatively 
high orbiting satellites as 1959 a1 and 1959 n, the effect - 
of the air drag cannot be described with sufficient 
accuracy in too long an interval of time by simple 
polynomial expressions. As a compromise, some nu- 
merical experimentation led to the 20-day intervals 
noted above, with third-degree polynomials in the 
mean anomaly M. The best possible polynomial repre- 
sentation of all the orbital elements was determined, 
treating w and Q as polynomials of the second degree, 
I and e as polynomials of the first degree, and M asa 
polynomial of the third degree. This means that the 
long-range direct and indirect effects of the air drag, 
of the radiation pressure, of secular and long-periodic 
perturbations due to the even and odd zonal harmonics 
in the earth’s gravitational potential, and of luni-solar 
perturbations were treated in an empirical way. Ana- 
lytical expressions for the first-order short-periodic 
perturbations due to the oblateness of the earth are 
incorporated in the computer program. The coefficients 
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TaBLeE II. Results of computations. 


Satellite 
1959 al 1959 » 
to April 16.0, 1960 May 19.0, 1960 
en 20421383 236°6646 
0; 1.002738 rev/day 
Q 138°0499 155°1312 
ont —3°51067/day —3°27567/day 
n 11.463937 rev/day 11.069014 rev/day 
a 1.302343 a, 1.333153 a, 
e 0.1646379 0.1891083 
q=a(i—e)—a, 561 km 517 km 
hi 32°8796 3323569 


of these polynomials were improved by including the 
trigonometric terms (8) into the orbital elements with 
the constants 


Bese <> anid A= O- 


‘of Heiskanen. In several consecutive runs these 
constants were then treated as unknowns only in the 
right ascension of the node , because one would 
expect no short-periodic air drag effects in this orbital 
element, and because the amplitude of 6Q turned out 
to be the biggest one. The improved values of 6 and \x 
were substituted into 6w, 62, 62, and 6M, and then the 
computation of these constants was repeated until the 
new values agreed with the old ones within their 
standard errors. The relevant results of these compu- 
tations are given in Table II. Consequently, expressions 
(8) become 


bw —BX0.4641 sin2 (@x—Q) 
62 8X 1.4808 sin2 (@x—) 
él 8X0.9572 cos2 (@x—Q) 
6M 8X0.7688 sin2 (@x—Q) 


—8X0.4041 sin2 (@x—2) 
8X1.3824 sin2 (6x—) 
8X0.9100 cos2 (@x—2) 
8X0.7370 sin2 (@x—Q). 


In the course of the trigonometric least-squares fit to 
the observations, the IBM 704 electronic computer 
actually obtained the amplitude of 6Q and twice the 
phase angle @x0—20=Ax+ (6g0—Q a From these data 
we immediately obtain 
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B (3.47-40.28) x 10-5 
Ax —33°5+1°9 


(2.88-+0.32) x 10-5 
Sey cea 


The agreement between the two sets of data seems 
reasonable; taking the weighted mean values, we obtain 


B= (3.210.29) 10-5, 
and 
Ax = —33°15+0°53. 


This numerical value of the ellipticity gives the per- 
turbation 6w, --- expressed in radians. For other angular 
units we note that 


B=5.11X10-* rev= 1°84 10-*=6762. 


According to the definition of the quantity 6, the differ- 
ence between the axes a, and 6, is 205 m. 

Needless to say, this determination of the coefficients 
of the second sectorial harmonic is far from being a 
final one. It will be improved as soon as more accurate 
coordinates of the Baker-Nunn camera stations are 
available. Also, the effect of other tesseral harmonics on 
this determination should be examined, possibly by 
using satellites with essentially different orbital elements. 
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Departures from the Russell Model in TX Ursae Majoris 
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Yellow, blue, and violet photoelectric light curves have been obtained for TX UMa. Additional obser- 
vations place the variable on the U, B, V system. It is demonstrated that the light curve is variable, at 
least at certain phases. The outside-eclipse variations exhibit departures from the Russell model as a 
function of wavelength of observation. Consistency in the solution of primary eclipse can be obtained if 
one abandons the Russell model for an ordinary body eclipse in favor of one comprising two sharp-edged 
stars, one of which is surrounded by an extensive but sharp-edged atmosphere. 
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I, INTRODUCTION 


HE most recent photometric studies of TX UMa 
are those by Wood (1946) and Huffer and Eggen 
(1947). Each of these articles adequately summarizes 
the past history of the eclipsing pair. Plavec (1960) has 
collected the published times of minimum light and 
studied the apsidal motion of the orbit. In addition, 
Struve and Huang (1956) have considered the possi- 
bility that the hot component is at times surrounded 
by an emitting ring composed of material ejected from 
the cool star. 
The notation used in the present paper for charac- 
teristics explicable by the Russell model is that recom- 
mended by Russell and Merrill (1952a). 


Il. OBSERVATIONS 


In 1956 and 1957 TX UMa was observed with the 
photoelectric photometers on the 36-inch and 12.5-inch 
reflectors of the Steward Observatory. A total of 372 
yellow, 365 blue, and 368 ultraviolet (hereafter, violet) 
observations were obtained. Those observations taken 
with the smaller reflector have been transformed to the 
instrumental system of the 36-inch telescope by the 
technique described by Koch (1960). All observations 
have been corrected for differential extinction using the 
mean coefficients given in the article just cited. The 
comparison star is BD +47°1797, and it was checked 
at least once in each color every night against BD 


(V-0), 
-0.5 
por sea, 
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¢ 
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Fic. 1. Yellow observations of TX UMa collected 
onto one light curve. 


+47°1799. The check star is ADS 7886, seemingly an 

optical double, and both components have been ob- | 
served through the aperture of the diaphragm slide. } 
Study of the magnitude differences between the com- } 
parison and check stars indicates that there are no | 
measurable night or seasonal variations in either star. — 
The observations are listed in Table I (a, b, and c) in | 
which time is counted from Greenwich Mean Noon. | 
The yellow observations are collected into one light | 
curve in Fig. 1. 


transformed to the U,B,V system of Johnson and © 
Morgan (1953) by using 6 UMa and Groombridge 1618 
as standards. The magnitudes and color indices for the 
comparison stars, determined by least-squares solu- 
tions, are shown in Table II. This table also lists the 
probable error for each determination and the number 
of observations of these stars. From these observations 
it is possible to scale the light variation of TX UMa on 
the U,B,V system, and this is shown in the last four 
lines of Table II. 


II. EPOCH AND PERIOD 


Newly determined -heliocentric times of minimum 
light in primary eclipse are JD 2435541.8634 and JD 
2435584.7500. These minima, together with all others’ 
published and available, have been tested against the 
ephemeris given by Wood. It is concluded that his 
light elements can be slightly improved to: 


Pr. Min. = 2416426.781-+-3.0633175 E+- 
+0.021 sin (0°09474 E+ 64°42). 


With these elements primary eclipse is symmetric. Since 
the orbital eccentricity is small, this is expected. 
Plavec’s ephemeris differs slightly from this one and 
shifts primary minimum from zero phase. Accordingly, 
the phases of the observations in Table I have been 
computed with the elements given above. The exact 
behavior of secondary eclipse in the new light curves 
is dealt with in a subsequent section. 


IV. VARIATIONS OUTSIDE ECLIPSES 


Night-by-night plots of the observations demon- 
strated that there are no consistent seasonal differences 
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TEX WEN: on 
TABLE Ia. Yellow observations of TX UMa. 
JD (hel.) - JD (hel.) JD (hel.) 
V-C 2435000-+ Phase V-C 2435000+ Phase V-C 2435000-+ Phase 
—0™351 528.8459 0.7507 —0™273 547.7884 0.9344 —0™283 555.6311 0.4945 
—0.345 .8581 0.7547 +0.049 .8561 0.9565 —0.272 . 6384 0.4969 
—0.331 . 8667 0.7575 +0.121 -8651 0.9594 —0.259 . 6488 0.5003 
—0.337 .8754 0.7604 +0.191 . 8728 0.9619 —0.275 .6557 0.5026 
—0.340 . 8845 0.7633 +0.325 . 8835 0.9654 —0.273 .6651 0.5056 
—0.355 530.6851 0.3511 +0.392 .8912 0.9679 —0.260 .6721 0.5079 
—0.324 . 6987 0.3556 +0.495 - 9009 0.9711 —0.274 .6814 0.5110 
—0.334 -7105 0.3594 +0.623 .9085 0.9736 —0.264 . 6894 0.5136 
—0.325 .7198 0.3625 +0.745 .9193 0.9771 —0.282 .6998 0.5170 
—0.329 .7285 0.3653 +0.868 -9269 0.9796 —0.283 .7071 0.5194 
—0.331 . 7372 0.3681 +1.034 .9380 0.9832 —0.286 .7179 0.5229 
—0.335 . 7459 0.3710 +1.123 .9467 0.9860 —0.277 .7262 0.5256 
—0.334 . 7556 0.3742 +1.273 -9582 0.9898 —0.288 .7370 0.5291 
—0.325 . 7684 0.3783 +1.348 . 9679 0.9930 —0.301 .7519 0.5340 
—0.333 nih 0.3812 +1.395 .9801 0.9969 —0.311 .7596 0.5365 
—0.331 . 7865 0.3842 +1.447 .9901 0.0002 —0.299 . 7693 0.5397 
—0.337 .7955 0.3872 —0.331 551.6677 0.1945 —0.311 .7790 0.5428 
—0.336 - 8063 0.3907 —0.341 . 6861 0.2067 —0.316 . 7887 0.5460 
—0.319 -8150 0.3935 —0.310 .6959 0.2099 —0.314. .7978 0.5490 
—(0.345 . 8247 0.3967 —0.326 . 7084 0.2140 —0.344 . 8096 0.5528 
—0.345 .8337 0.3996 —0.347 .7170 0.2168 —0.333 .8169 0.5552 
—0.343 . 8441 0.4030 —0.343 . 7282 0.2205 —0.326 .8273 0.5586 
—0.343 .8577 0.4075 —0.326 . 7354 0.2228 —0.327 . 8346 0.5613 
—0.336 . 8674 0.4106 —0.320 . 7462 0.2263 —0.345 8443 0.5641 
—0.333 .8778 0.4140 —0.340 -7743 0.2355 —0.339 .8519 0.5666 
—0.333 .8917 0.4186 —0.329 . 7809 0.2377 —0.333 . 8620 0.5699 
—0.336 -9014 0.4217 —0.335 . 7900 0.2406 —0.335 . 8693 0.5723 
—0.341 -9136 0.4257 —0.349 . 8386 0.2565 —0.339 .8790 0.5755 
—0.334 .9278 0.4304 —0.354 .8479 0.2595 —0.341 - 8863 0.5779 
—0.340 .9389 0.4340 —0.345 .8552 0.2619 —0.339 . 8974 0.5815 
—0.330 .9490 0.4373 —0.279 552.6201 0.5116 —0.311 556.7164 0.8488 
—0.326 -9584 0.4404 —0.271 .6347 0.5164 —0.313 .7310 0.8536 
—0.327 -9702 0.4442 —0.269 .6417 0.5187 —0.317 - 7383 0.8560 
+1.187 538.8366 0.0121 —0.274 .6510 0.5217 —0.319 7477 0.8591 
+1.040 . 8495 0.0163 —0.292 -6597 0.5245 —0.322 . 7550 0.8614 
+0.880 .8613 0.0202 —0.286 . 6698 0.5278 —0.313 . 7664 0.8652 
+0.766 - 8696 0.0229 —0.284 .6771 0.5302 —0.307 . 7748 0.8679 
+0.576 - 8839 0.0276 —0.297 .6861 0.5332 —0.316 - 7855 0.8714 
+0.444 . 8953 0.0313 —0.304 .6927 0.5353 —0.322 . 7946 0.8744 
+0.299 - 9068 0.0350 —0.310 .7021 0.5384 —0.319 . 8039 0.8774 
+0.010 -9179 0.0387 —0.298 . 7090 0.5406 —0.315 .8116 0.8799 
—0.023 -9412 0.0463 —0.306 .7156 0.5428 —0.308 .8241 0.8840 
—0.093 -9523 0.0499 —0.318 (GAME 0.5461 —0.315 . 8328 0.8868 
—0.195 -9738 0.0569 —0.322 .7385 0.5503 —0.320 557.7720 0.1934 
—0.186 541.6931 0.9446 —0.320 . 7496 0.5539 —0.322 . 7858 0.1979 
+0.152 .7425 0.9607 —0.319 -7566 0.5562 —0.323 .7928 0.2002 
+0.215 .7491 0.9629 —0.322 . 7667 0.5595 —0.326 .8015 0.2031 
+0.334 .7591 0.9662 —0.324 .7739 0.5618 —0.318 .8091 0.2055 
+0.456 .7702 0.9698 —0.318 . 7837 0.5650 —0.325 559.7541 0.8405 
+0.626 . 7831 0.9740 —0.340 -7910 0.5674 —0.310 . 7684 0.8451 
+1.088 .8181 0.9854 —0.326 . 8007 0.5706 —0.315 ise) 0.8474 
+1.255 - 8293 0.9891 —0.324 .8076 0.5728 —0.326 . 7854 0.8507 
+1.364 .8376 0.9918 —0.289 554.6200 0.1645 —0.327 . 7934 0.8533 
+1.392 . 8497 0.9957 —0.312 - 6329 0.1689 —0.313 .8027 0.8563 
+1.435 . 8584 0.9986 —0.310 .6395 0.1708 —0.309 . 8093 0.8585 
+1.419 .8727 0.0032 —0.321 . 6488 0.1739 —0.257 569.6205 0.0612 
+1.386 . 8810 0.0059 —0.325 .6544 0.1757 —0.295 -6365 0.0665 
+1.272 .8931 0.0099 —0.320 .6655 0.1793 —0.294 .6435 0.0688 
+1.192 -9015 0.0126 —0.320 .6721 0.1815 —0.291 -6525 0.0717 
+1.043 -9133 0.0165 —0.326 . 6839 0.1853 —0.300 .6612 0.0745 
+0.924 -9220 0.0193 —0.325 . 6902 0.1874 —0.287 -6705 0.0776 
+0.745 -9341 0.0233 —0.327 .6992 0.1903 —0.283 .6771 0.0797 
+0.605 -9435 0.0263 —0.323 .7058 0.1925 —0.302 . 7063 0.0893 
—0.349 542.8145 0.3107 —0.325 .7159 0.1958 —0.302 SPELL 0.0928 
—0.350 .8244 0.3139 —0.334 .7291 0.2001 —0.299 .7240 0.0950 
—0.356 543.8844 0.6599 —0.333 . 7388 0.2033 —0.292 . 7334 0.0981 
—0.341 .8972 0.6641 —0.329 . 7495 0.2068 —0.297 - 7417 0.1008 
—0.338 - 9083 0.6677 —0.322 .7565 0.2090 —0.300 STV 0.1042 
—0.340 -9160 0.6702 —0.337 . 7662 0.2122 —0.294 . 7876 0.1158 
—0.348 .9271 0.6739 —0.334 AG eH: 0.2145 —0.299 .7966 . 0.1187 
—0.337 -9347 0.6763 —0.335 . 7846 0.2182 —0.338 579.6431 0.3331 
—0.347 - 9448 0.6796 —0.337 -7923 0.2207 —0.327 .6588 0.3382 
—0.351 -9552 0.6830 —0.276 555.6175 0.4901 —0.329 . 6865 0.3473 
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JD (hel.) JD (hel.) JD (hel.) 
vV-C 2435000-+- Phase Vee 2435000-++ Phase V-C€ 2435000-++ 
—0™342 579.6970 0.3507 =02331 818.9257 0.4449 —0™316 860.8849 
+0.336 584. 6440 0.9656 —0.326 -9330 0.4472 —0.339 862.7871 
+0.550 6645 0.9723 —0.319 .9424 0.4503 —0.335 7957 
+0.728 6782 0.9767 —0.314 9493 0.4526 —0.335 - 8044 
+0.841 . 6865 0.9794 —0.306 9663 0.4581 —0.325 -8131 
+1.002 6969 0.9828 —0.313 9730 0.4603 —0.318 8218 
+1.108 . 7059 0.9858 —0.309 .9830 0.4636 —0.306 893.8262 
+1.240 7178 0.9897 —0.296 9903 0.4659 —0.286 - 8404 
S15 289 7254 0.9921 —0.302 819.0010 0.4694 —0.291 8491 
+1.354 . 7368 0.9958 —0.286 .0104 0.4725 —0.291 8578 
+1.389 7441 0.9982 —0.291 .0278 0.4782 —0. 288 . 8668 
+1.389 . 7563 0.0022 —0.282 .0361 0.4809 —0.293 -8758 
+1.340 . 7671 0.0058 —0.334 8931 0.7607 0.284 - 8848 
= ESP sis) .7778 0.0093 —0.330 9035 0.7644 —0.276 .9057 
+1.207 . 7848 0.0115 —0.341 .9119 0.7668 —0.289 .9157 
+1.052 . 7955 0.0150 —0.309 .9296 0.7726 —0.288 9251 
+0.961 8039 0.0178 —0.321 .9362 0.7747 —0.287 9345 
+0.698 8223 0.0238 —0.328 .9456 0.7778 —0.289 -9446 
+0.564 .8316 0.0268 —0.326 -9522 0.7800 —0.268 .9533 
+0.452 8428 0.0305 —0.328 .9619 0.7831 —0.238 .9626 
+0.339 8507 0.0331 —0.326 9685 0.7853 —0.199 .9713 
+0.155 - 8688 0.0390 —0.322 .9778 0.7883 —0.157 9803 
+0.066 8771 0.0417 —0.325 9844 0.7905 —0.104 9893 
+0.003 -8879 0.0452 —0.320 .9942 0.7937 —0.037 .9991 
—0.334 588.7685 0.3120 —0.295 820.9161 0.0946 +0.018 894.0084 
—0.336 7831 0.3167 —0.312 “9484 0.1052 —0.336 931.7218 
—0.343 7911 0.3193 —0.290 -9595 0.1088 —0.335 SSO 
—0.321 8025 0.3231 —0.308 9672 0.1113 —0.342 . 7444 
—0.341 -8105 0.3257 —0.315 .9793 0.1153 —0.329 -7530 
—0.330 8213 0.3292 —0.321 9925 0.1196 —0.330 .7617 
—0.369 8313 0.3325 —0.302 821.0026 0.1229 —0.343 . 7704 
—0.355 589.6695 0.6061 —0.336 837.8065 0.6083 —0.312 946. 6444 
—0.290 615.6486 0.0868 =0:319 840.8785 0.6112 —0.316 6604 
—0.296 6643 0.0919 —0.330 -8875 0.6141 —0.309 6698 
—0.290 6733 0.0948 —0.327 8945 0.6164 —0.306 6784 
—0.346 617.6554 0.7419 —0.330 -9035 0.6193 —0.289 6874 
—0.328 -6714 0:7471 —0.330 9098 0.6214 —0.301 6965 
=0.325 .6790 0.7496 —0.327 -9184 0.6242 —0.342 954.7295 
—0.328 6887 0.7528 —0.328 -9250 0.6263 —0.339 . 7420 
—0.349 6981 0.7558 —0.328 9341 0.6293 —0.345 .7514 
—0.336 7342 0.7676 —0.330 .9410 0.6316 —0.333 .7611 
—0.325 - 7443 0.7709 —0.331 -9695 0.6409 —0.336 . 7702 
Ooo HOS 0.7740 =0.319 .9792 0.6440 —0.339 7795 
+0.626 624.6535 0.0263 —0.338 .9862 0.6463 —0.330 . 7886 
+0. 388 -6719 0.0323 —0.346 9959 0.6495 . —0.335 71973 
+0. 290 6823 0.0357 —0.349 841.0032 0.6519 —0.331 -8059 
+0.191 6931 0.0393 =0-337 .0132 0.6551 —0.331 .8164 
+0.063 - 1028 0.0424 —0.301 860.8446 0.1289 —0.315 8250 
O29 818.8958 0.4351 OST .8561 OMS 7 —0.335 8351 
—0.329 9097 0.4396 —0.326 -8655 0.1358 —0.332 8441 
—0.334 9163 0.4418 —0.304 8731 0.1382 —0.330 8532 


within the light curves. On JD 2435818-9 and JD 
2435893-4, however, the observations are respectively 
systematically high and low in a suspicious manner. 
Nonetheless, it was possible to analyze the observations 
outside eclipse for the expected Fourier expansions 
according to the graphical method of Russell and 
Merrill (1952b). These are given as 


£(y)=0.9735—0.0192 cos#—0.0139 cos26, 
£(b)=0.9819—0.0139 cos6—0.0056 cos26, 
£(u)=0.9780—0.0162 cosd—0.0058 cos26. 


The probable residuals for a single observation from 
these computed variations are +0™006 in yellow and 


blue and +0008 in violet. Terms in cos@ higher than 
those listed and all terms in sin@ are small enough to be 
safely neglected. 


V. PRELIMINARY RECTIFICATIONS AND SOLUTIONS 


From the solutions given previously by O’Keefe 
(1938) and Wood it is possible to derive approximate 
values of the inclination of the orbit to the plane of the 
sky and the mean ratio of the surface brightnesses of 
the components as a function of wavelength. With these 
quantities known, the observations were provisionally 
rectified for reflection and the mean ellipticity of figure 
of the stars according to the following equations: 


TABLE Ib. Blue observations of TX UMa. 


se, Gn O00 


JD (hel.) JD (hel.) JD (hel.) 
V-C 2435000-+ Phase V-C 2435000+ Phase V-C 2435000+- Phase 
—0"514 528.8480 0.7514 +0007 547.8658 0.9596 —0™430 555.6478 0.5000 
—0.508 . 8574 0.7545 +0.095 .8735 0.9621 —0.451 -6547 0.5022 
—0.494 . 8660 0.7573 +0.229 . 8842 0.9656 —0.442 .6641 0.5053 
—0.487 .8747 0.7601 +0.318 .8919 0.9681 —0.446 .6710 0.5076 
—0.497 - 8838 0.7631 +0.459 .9016 0.9713 —0.441 . 6804 0.5106 
—0.514 530.6844 0.3509 +0.591 .9092 0.9738 —0.440 . 6884 0.5133 
—0.484 . 6980 0.3554 +0.766 .9200 0.9773 —0.456 -6988 0.5166 
—0.488 . 7098 0.3592 +0.909 -9276 0.9798 —0.439 .7061 0.5190 
—0.486 .7191 0.3622 +1.117 - 9387 0.9834 —0.451 .7169 0.5226 
—0.494 .7278 0.3651 +1.249 .9474 0.9863 —0.450 252 0.5253 
—0.497 . 7365 0.3679 +1.429 9589 0.9900 —0.453 .7360 0.5288 
—0.498 .74140 0.3694 +1.557 .9686 0.9932 —0.459 . 7509 0.5337 
—0.495 . 7549 0.3739 +1.640 - 9808 0.9972 —0.473 .7585 0.5361 
—0.489 . 7678 0.3781 +1.661 .9908 0.0004 —0.468 . 7682 0.5393 
—0.483 .7764 0.3809 —0.470 551.6670 0.2005 —0.473 .7780 0.5425 
—0.490 . 7858 0.3840 —0.495 - 6854 0.2065 —0.472 -1877 0.5457 
—0.489 . 7948 0.3869 —0.461 .6952 0.2097 —0.466 . 7946 0.5479 
—0.498 .8056 0.3905 —0.491 . 7098 0.2145 —0.500 .8085 0.5526 
—9.489 .8143 0.3933 —0.490 . 7164 0.2166 —0.485. .8158 0.5548 
—0.500 . 8240 0.3965 —0.483 :1215 0.2202 —0.470 . 8262 0.5582 
—0.500 - 8330 0.3994 —0.485 . 7348 0.2226 —0.478 - 8335 0.5606 
—0.514 . 8434 0.4028 —0.483 . 1455 0.2261 —0.479 . 8432 0.5638 
—0.510 .8570 0.4073 —0.481 . 7736 0.2353 —0.492 . 8509 0.5663 
—0.483 . 8667 0.4104 —0.480 . 7802 0.2374 —0.496 . 8606 0.5695 
—0.488 8771 0.4138 —0.476 - 7893 0.2404 —0.501 . 8682 0.5719 
—0.473 .8875 0.4172 —0.494 . 8379 0.2563 —0.493 .8780 0.5751 
—0.493 - 9007 0.4215 —0.490 - 8473 0.2593 —0.498 -8853 0.5775 
—0.494 -9129 0.4255 —0.487 .8545 0.2617 —0.483 . 8964 0.5812 
—0.498 .9271 0.4301 —0.456 552.6187 0.5112 —0.476 556.7154 0.8485 
—0.498 .9382 | 0.4338 —0.461 . 6337 0.5161 —0.477 .7300 0.8533 
—0.480 . 9483 0.4371 —0.457 . 6406 0.5183 —0.473 .1373 0.8557 
—0.485 .9577 0.4401 —0.452 . 6500 0.5214 —0.470 . 7467 0.8587 
—0.492 -9695 0.4440 —0.464 6618 0.5252 —0.478 .7539 0.8611 
+1.326 538.8373 0.0124 —0.461 . 6687 0.5275 —0.472 - 7633 0.8641 
+1.112 . 8502 0.0166 —0.456 . 6760 0.5299 —0.464 .7730 0.8673 
+0.887 - 8620 0.0204 —0.467 .6851 0.5328 —0.476 . 7845 0.8711 
+0.755 -8703 0.0231 —0.475 .6917 0.5350 —0.473 -7935 0.8740 
+0.517 . 8846 0.0278 —0.467 . 7010 0.5380 —0.478 . 8029 0.8771 
+0.356 . 8960 0.0315 —0.468 . 7080 0.5403 —0.475 .8105 0.8796 
+0.185 -9075 0.0353 —0.471 . 7146 0.5425 —0.474 . 8230 0.8836 
+0.003 -9186 0.0389 —0.472 . 7246 0.5457 —0.477 .8317 0.8865 
—0.179 .9418 0.0465 —0.492 RAS 0.5499 —0.470 557.7695 0.1926 
—0.235 9516 0.0497 —0.480 . 7486 0.5536 —0.474 . 7848 0.1976 
—0.354 .9731 0.0567 —0.484 WEES 0.5558 —0.477 7917 0.1999 
+0.025 541.7418 0.9605 —0.480 .7656 0.5591 —0.471 . 8004 0.2027 
+0.123 7497 0.9631 —0.480 .7729 0.5615 —0.468 . 8081 0.2052 
+0.245 . 7598 0.9664 —0.471 . 7826 0.5647 —0.478 559.7531 0.8402 
+0.402 .7709 0.9700 —0.481 . 7899 0.5670 —0.472 . 1673 0.8448 
+0.609 . 7838 0.9742 —0.474 . 7996 0.5702 —0.481 .1743 0.8471 
+1.230 . 8188 0.9856 —0.485 .8066 0.5725 —0.477 . 7843 0.8503 
+1.409 . 8300 0.9893 —0.467 554.6190 0.1642 —0.483 . 7923 0.8529 
+1.510 .8383 0.9920 —0.463 .6318 0.1683 —0.477 .8017 0.8560 
+1.643 - 8504 0.9960 —0.472 .6384 0.1705 —0.480 . 8083 0.8582 
+1.686 . 8591 0.9988 —0.472 .6478 0.1736 —0.472 569.6355 0.0662 
+1.651 .8734 0.0034 —0.477 .6544 0.1757 —0.467 . 6424 0.0684 
+1.609 . 8817 0.0062 —0.480 . 6645 0.1790 —0.453 .6515 0.0714 
+1.445 . 8938 0.0101 —0.480 .6711 0.1812 —0.477 .7053 0.0889 
+1.314 .9022 0.0128 —0.479 . 6801 0.1841 —0.465 . 7160 0.0924 
+1.122 -9140 0.0167 —0.485 . 6891 0.1870 —0.464 . 7230 0.0947 
+0.957 -9227 0.0195 —0.482 . 6982 0.1900 —0.463 . 7324 0.0978 
+0.728 -9348 0.0235 —0.483 . 7048 0.1922 —0.459 . 7407 0.1004 
+0.551 -9442 0.0266 —0.488 .7148 0.1954 —0.465 1511 0.1039 
—0.488 542.8138 0.3104 —0.476 .7280 0.1997 —0.461 . 7865 0.1144 
—0.505 _ .8237 0.3137 —0.476 RT! 0.2029 —0.474 . 7938 0.1178 
—0.498 543.8837 0.6597 —0.470 . 7485 0.2064 —0.482 579.6421 0.3328 
—0.489 .8965 0.6639 —0.475 . 7554 0.2087 —0.489 .6577 0.3379 
—0.478 -9077 0.6675 —0.487 . 7652 0.2119 —0.493 . 6855 0.3469 
—0.479 -9153 0.6700 —0.487 AeA! 0.2141 —0.495 .6959 0.3503 
—0,498 -9264 0.6736 —0.481 . 7836 0.2179 +0.231 584.6443 0.9657 
—0.482 -9340 0.6761 —0.486 .7912 0.2204 +0.456 . 6613 0.9712 
—0.483 -9441 0.6794 —0.467 555.6165 0.4898 +0.764 -6813 0.9777 
—0.473 -9511 0.6817 —0.455 .6300 0.4942 +1.073 .6976 0.9831 
—0.065 547.8575 0.9569 —0.444 .6373 0.4966 +1,222 7066 0.9860 


234 RIO iB EAR Een OF CET 
TasE Ib.—Continued. 
JD (hel.) JD (hel.) JD (hel.) 
Vv-C 2435000-+- Phase V-C 2435000++ Phase VEC 2435000-++ 
+17406 584.7184 0.9899 —0™487 818.9737 0.4605 —0™494 862.8051 
+1.487 .7261 0.9924 —0.487 .9837 0.4638 —0.479 -8138 
+1.604 Sets 0.9961 —0.463 .9910 0.4662 —0.486 8225 
+1.639 . 7448 0.9985 —0.477 819.0017 0.4697 —0.459 893 . 8269 
+1.657 7591 0.0032 —0.474 -O111 0.4727 —0.449 .8411 
+1.560 . 7678 0.0060 —0.463 .0285 0.4784 —0.449 8498 
+1.449 7785 0.0095 —0.461 .0368 0.4812 —0.453 - 8584 
+1.343 . 7855 0.0118 —0.482 .8938 0.7609 —0.451 -8675 
+1.164 . 7962 0.0153 —0.475 .9042 0.7643 —0.445 .8765 
+1.028 . 8046 0.0180 —0.487 -9126 0.7670 —0.448 .8855 
+0.702 .8230 0.0240 —0.468 -9303 0.7728 —0.451 -9064 
+0.548 8323 0.0270 —0.473 .9369 0.7750 —0.445 -9164 
-+-0.370 . 8434 0.0307 —0.483 -9462 0.7780 ~—0.448 .9258 
+0.251 .8514 0.0333 —0.491 -9528 0.7802 —0.442 -9352 
+0.046 .8678 0.0385 =—0.480 .9626 0.7834 —0.436 .9453 
—0.042 .8760 0.0413 —0.478 -9692 0.7855 —0.415 -9539 
—0.127 .8868 0.0448 —0.477 -9785 0.7885 —0.399 -9633 
—0.492 588 . 7692 0.3122 —0.478 ~ -9851 0.7907 —0.363 .9720 
—0.492 . 7838 0.3170 —0.477 .9949 0.7939 —0.307 .9810 
—0.488 .7918 0.3196 —0.475 820.9168 0.0948 —0.252 .9900 
—0.480 . 8032 0.3233 —0.479 9491 0.1054 —0.195 .9998 
—0.496 .8112 0.3259 —0.467 .9602 0.1090 —0.119 894.0091 
—0).492 .8220 0.3294 —0.479 .9678 0.1115 —0.491 931.7225 
—0.523 8306 0.3322 —0.478 .9800 0.1155 —0.493 1364 
—0.502 589.6702 0.6063 —0.483 5 BS) 0.1199 —0.493 - 7451 
—0.443 615.6493 0.0870 —0.464 821.0033 0.1231 —0.482 (eit 
—0.455 . 6650 0.0921 —0.498 837.8072 0.6086 —0.489 .7624 
—0.475 .6740 0.0951 —0.483 840.8792 0.6114 —0.487 Ea 
—0.488 617.6547 0.7417 —0.486 - 8882 0.6143 —0.471 946.6451 
—0.482 .6721 0.7473 —0.481 -8952 0.6166 —0.479 .6611 
—0.464 .6797 0.7498 —0.486 -9042 0.6195 —0.467 -6704 
—0.475 .6908 0.7534 —0.488 .9105 0.6216 —0.468 .6791 
—0.491 6988 0.7561 —(.481 .9191 0.6244 —0.438 6381 
—0.481 . 7349 0.7678 —0.479 -9257 0.6266 —0.456 -6972 
—0.491 . 7450 0.7711 —0.491 -9348 0.6295 —0.491 954.7302 
—0.492 . 7544 0.7742 —0.491 -9417 0.6318 —0.490 7427 
+0.540 624.6559 0:0271 —0.487 .9702 0.6411 —0.491 7521 
+0.330 .6712 0.0321 —0.475 .9799 0.6443 —0.492 .7618 
+0.196 -6813 0.0354 —0.487 -9869 0.6465 —0.479 .7709 
+0.069 .6920 0.0389 —0.495 -9966 0.6497 —0.483 . 7802 
—0.010 . 7007 0.0417 —0.498 841.0039 0.6521 —0.495 . 7893 
—0.473 818.8965 0.4353 —0.493 .0139 0.6554 —0.490 .7979 
—0.495 -9104 0.4399 —0.476 860.8453 0.1292 —0.486 - 8066 
—0.494 .9170 0.4420 —0.459 .8547 0.1322 —0.484 .8170 
—0.488 .9264 0.4451 —0.478 8641 0.1353 -—0.486 -8257 
—0.480 .9337 0.4475 —0.472 .8738 0.1385 —0.483 -8358 
—0.475 -9431 0.4505 —0.479 - 8856 0.1423 —0.460 - 8448 
—0.470 -9500 0.4528 —0.485 862.7878 0.7633 —0.503 - 8538 
—0.464 .9670 0.4583 —0.487 * 7964 0.7661 
s’G)= £(y)+0.0156-+0.0192 cosd+-0.0051 cos20 were a circle, but as has been proved by Russell (1912), 
y 0.9891 —0.0088 cos26 , the departures from them, if strict account be taken of 
the eccentric orbit, are very small. 
N2z(y)=0.0778, It is well known that, except for the particular orien- 
j tation with the major axis in the plane of the sky, the 
g! ee (6) +-0.0120+-0.0139 cos#+-0.0039 cos26 exact and complete solution of a light curve for stars 
0.9939—0.0017 cos26 2 moving in an eccentric orbit proceeds by successive 
approximations. The details of the method were given 
N2z(b)=0.0388, first by Russell and subsequently restated and some- 
what elaborated by Kopal (1946). The technique 
9” (u)= &(u)+0.0153-+0.0162 cos6-+0.0050 cos26 involves the pairing up of the geometric and scale 


) 


0.9933—0.0008 cos26 
Nz(u)=0.0444. 


These equations are derived as if the orbit of TX UMa 


elements which can be derived individually from 
primary and sufficiently deep secondary minima so as 
to obtain a consistent solution. To the degree of accur- 
acy that quantities of the second order in the orbital 
eccentricity can be neglected, the equations for con- 


TX UMa 235 
TABLE Ic. Violet observations of TX UMa. 
JD (hel.) JD (hel.) JD (hel.) 
WC 2435000-+- Phase vV—€ 2435000+- Phase WAC 2435000+- Phase 
—0™"945 528.8466 0.7510 — 0394 547.8665 0.9599 —O™885 555.6495 0.5006 
—0.949 . 8588 0.7549 —0.285 .8742 0.9624 —0.899 .6564 0.5028 
—0.935 .8674 0.7578 —0.144 . 8849 0.9659 —0.903 .6658 0.5059 
—0.932 .8761 0.7606 —0.054 .8926 0.9684 —0.890 .6728 0.5082 
—0.945 .8851 0.7635 +0. 106 -9023 0.9715 —0.914 -6821 0.5112 
—0.945 530.6858 0.3514 +0.267 .9099 0.9740 —0.897 .6901 0.5138 
—0.906 .6994 0.3558 +0.477 -9207 0.9776 —0.916 . 7005 0.5172 
—0.923 .7112 0.3597 +0.638 -9283 0.9800 —0.902 -7078 0.5196 
—0.921 .7205 0.3627 +0.885 .9394 0.9837 —0.892 .7200 0.5236 
—0.934 .7292 0.3655 +1.019 .9481 0.9865 —0.905 .7269 0.5258 
—0.936 .7379 0.3684 +1.274 -9596 0.9903 —0.911 BOM 0.5293 
—0.928 . 7466 0.3712 +1.453 .9693 0.9934 —0.931 .7526 0.5342 
—0.937 . 7563 0.3744 +1.565 -9814 0.9974 —0.941 .7603 0.5367 
—0.927 -7691 0.3786 +1.592 -9915 0.0007 —0.925 .7700 0.5399 
—0.918 .7778 0.3814 —0.923 551.6684 0.2009 —0.933 .7797 0.5431 
—0.922 . 7872 0.3845 —0.941 . 6868 0.2069 —0.947 . 7894 0.5462 
—0.930 . 7962 0.3874 —0.900 .6966 0.2101 —0.934 . 7984 0.5492 
—0.942 .8070 0.3909 —0.906 .7091 0.2142 —0.945 .8103 0.5530 
—9.937 .8157 0.3938 —0.930 .7177 0.2170 —0.938 . .8175 0.5554 
—0.956 .8254 0.3969 —0.938 . 7288 0.2207 —0.934 .8280 0.5588 
—0.943 .8344 0.3999 —0.937 .7361 0.2230 —0.931 8353 0.5612 
—0.933 . 8448 0.4033 —0.928 . 7469 0.2266 —0.937 .8450 0.5644 
—0.946 . 8584 0.4077 —0.935 .7750 0.2357 —0.942 .8526 0.5669 
—0.948 . 8681 0.4109 —0.923 .7816 0.2379 —0.943 .8627 0.5702 
—0.930 .8785 0.4143 —0.930 . 7907 0.2409 —0.945 . 8700 0.5725 
—0.936 .8924 0.4188 —0.931 .8393 0.2567 —0.938 .8797 0.5757 
—0.935 -9021 0.4220 —0.942 . 8486 0.2598 —0.940 .8870 0.5781 
—0.951 -9143 0.4260 —0.945 .8559 0.2622 —0.941 .8981 0.5817 
—0.951 -9285 0.4306 —0.894 552.6194 0.5114 —0.905 556.7171 0.8491 
—0.934 .9396 0.4342 —0.911 .6354 0.5166 —0.936 .7317 0.8538 
—0.931 -9497 0.4375 —0.903 .6424 0.5189 —0.919 .7390 0.8562 
—0.927 .9591 0.4406 —0.889 .6517 0.5219 —0.921 . 7484 0.8593 
—0.945 .9723 0.4449 —0.925 . 6604 0.5248 —0.929 .7557 0.8617 
+1.109 538.8380 0.0126 —0.914 .6705 0.5281 —0.913 .7671 0.8654 
+0. 833 .8509 0.0168 —0.914 .6778 0.5305 —0.908 .7755 0.8681 
+0.581 . 8627 0.0206 —0.917 . 6868 0.5334 —0.920 . 7862 0.8716 
+0.413 .8710 0.0234 —0.929 .6934 0.5355 —0.910 .7953 0.8746 
+0.147 .8853 0.0280 —0.928 .7028 0.5386 —0.923 .8046 0.8776 
—0.047 . 8967 0.0317 —0.925 . 7097 0.5409 —0.914 .8102 0.8795 
—0.216 -9082 0.0355 —0.922 .7163 0.5430 —0.912 .8248 0.8842 
—0.267 -9193 0.0391 —0.931 .7264 0.5463 —0.924 8335 0.8871 
—0.608 -9425 0.0467 —0.936 .7392 0.5505 —0.900 557.7726 0.1936 
—0.677 -9530 0.0501 —0.936 . 7503 0.5541 —0.909 - 7865 0.1982 
—0.805 -9745 0.0571 —0.935 S7/S¥K8 0.5564 —0.907 .7935 0.2006 
—0.779 541.6938 0.9448 —0.931 .7674 0.5597 —0.931 .8022 0.2033 
—0.366 .7431 0.9609 —0.928 . 7746 0.5621 —0.906 559.7548 0.8407 
—0.264 . 7504 0.9633 —0.912 . 7844 0.5653 —0.907 .7691 0.8454 
—0.123 . 7605 0.9666 —0.940 .7917 0.5676 —0.908 .7760 0.8476 
+0.039 .7716 0.9702 —().929 .8014 0.5708 —0.928 .7861 0.8509 
+0.274 . 7845 0.9745 —0.930 . 8083 0.5731 —0.922 .7941 0.8536 
+1.018 .8195 ~ 0.9859 —0.910 554.6207 0.1647 —0.916 . 8034 0.8566 
+1.254 . 8306 0.9895 —0.901 -6336 0.1689 —0.912 .8100 0.8587 
+1.385 .8390 0.9922 —0.895 -6402 0.1711 —0.876 569.6212 0.0614 
+1.534 .8511 0.9962 —0.924 -6495 0.1741 —0.904 .6372 0.0667 
+1.594 . 8598 0.9990 —0.922 -6561 0.1763 —0.900 -6442 0.0690 
+1.546 .8741 0.0037 —0.936 . 6662 0.1796 —0.890 -6532 0.0719 
+1.484 . 8824 0.0064 —0.926 .6728 0.1817 —0.894 .6619 0.0748 
+1.279 . 8945 0.0103 —0.925 - 6846 0.1856 —0.886 .6712 0.0778 
+1. 108 -9029 0.0131 —0.922 . 6909 0.1876 —0.880 .6778 0.0800 
+0.849 .9147 0.0169 —0.918 . 6999 0.1906 —0.903 .7070 0.0895 
+0.659 .9234 0.0198 —0.911 .7065 0.1927 —0.912 .7178 0.0930 
—0.949 542.8152 0.3109 —0.913 .7166 0.1960 —0.912 .7247 0.0953 
—0.947 .8251 0.3141 —0.937 .7298 0.2003 —0.898 .7341 0.0983 
—0.934 543.8851 0.6602 —0.903 .7395 0.2033 —0.895 .7424 0.1011 
—0.922 .8979 0.6643 —0.905 .7502 0.2070 —0.906 .7528 0.1044 
—0.905 .9090 0.6680 —0.924 Lok 0.2093 —0.900 . 7883 0.1160 
—0.910 .9167 0.6705 —0.923 . 7669 0.2124 —0.912 .7973 0.1190 
—0.929 .9278 0.6741 —0.930 .7738 0.2147 —0.941 579.6449 0.3337 
—0,919 .9354 0.6766 —0.917 . 7853 0.2184 —0.940 -6595 0.3384 
—0.934 -9455 0.6799 —0.925 .7929 0.2209 —0.946 -6872 0.3475 
—0.921 .9545 0.6828 —0.909 555.6182 0.4903 —0.958 .6976 0.3509 
—0.898 547.7891 0.9346 —0.906 -6318 0.4948 —0.111 584.6471 0.9666 
—0.512 - 8568 0.9567 —0.914 -6391 0.4972 +0.149 -6641 0.9721 
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TABLE Ic.—Continued. 
JD (hel.) JD (hel.) 7 JD (hel.) 

VEC 2435000-+- Phase WES: 2435000-- Phase V-€ 2435000-- Phase 
+07455 584. 6803 0.9774 —0"938 818.9507 0.4530 — 02041 862.7971 0.7663 
+0.598 6879 0.9799 —0.936 9677 0.4586 —0.947 .8058 0.7691 
+0. 842 6983 0.9833 —0.941 9744. 0.4608 —0.927 8145 0.7720 
+1.091 . 7087 0.9867 —0.950 9844 0.4640 —0.935 8232 0.7748 
dedi . 7191 0.9901 —0.940 9917 0.4664 0.921 893.8275 0.8959 
+1.443 . 1282 0.9931 —0.947 819.0024 0.4699 —0.898 _ 8418 0.9006 . 
=f o30 . 7382 0.9963 —0.945 .0118 0.4729 —0.904 8505 0.9034. 
+1.559 7455 0.9987 —=05939 0292 0.4786 —0.907 8591 0.9062 
+1.609 SSV iT 0.0027 —0.936 0375 0.4814 —0.910 . 8682 0.9092 
+1.453 . 7684 0.0062 — 01925 8945 0.7611 —0.912 8772 0.9121 
+1.307 7792 0.0097 —0.928 9049 0.7645 —0.905 -8862 0.9151 
+1.156 . 7862 0.0120 —0.928 9133 0.7673 —0.906 9071 0.9219 
+0.914 7969 0.0155 0917 9310 0.7730 —0.902 9171 0.9252 
+0.741 8053 0.0182 =0; 927 9376 0.7752 ~—0.898 9265 0.9282 
+0.354 8237 0.0242 —0.925 9469 0.7782. —0.898 9359 0.9313 
+0.162 . 8330 0.0273 =0:933 9535 0.7804 —0.885 9459 0.9346 
—0.013 8441 0.0309 —0.937 9633 0.7836 —0.862 9546 0.9374 
—0.148 8521 0.0335 =0.933 -9699 0.7857 —0.827 -9640 0.9405 
—0.371 . 8695 0.0392 —0.930 19792 0.7888 —0.788 9727 0.9433 
—0.471 .8778 0.0419 =0.931 9858 0.7909 —0.729 9817 0.9462 
O99 . 8886 0.0454 =(0).935 9956 0.7941 —O.0a0 9907 0.9492 
—0.945 588.7699 0.3124 —0.898 820.9175 0.0951 —0.588 894.0005 0.9524 
—0.949 7845 0.3172 —0.927 9498 0.1056 —0.516 0098 0.9554 
—0.947 7925 0.3198 —0.893 9609 0.1092 —0.935 931.7232 0.2666 
—0.928 .8039 0.3235 =0°917 9685 0.1117 —0.941 hol 0.2712 
—0; 907, 8119 0.3261 = 5015 -9807 0.1157 —0.946 7457 0.2740 
=0:959 8237 0.3300 —0.926 * .9942 0.1201 —0.934 7544 0.2768 
—= 02070. 8320 0.3327 —0.923 821.0040 0.1236 —0.940 7631 0.2796 
—0.960 589.6709 0.6065 —0.938 837.8079 0. 6088 —0.943 1718 0.2825 
—0.903 615.6500 0.0872 =-O.029 840.8799 0.6116 —0.922 946.6458 0.1380 
—0.943 6656 0.0923 —0.934 . 8889 0.6145 —0.914 6618 0.1432 
—0.926 6761 0.0958 —0.944 8959 0.6168 —0.918 -6711 0.1462 
—0.950 617.6561 0.7421 —0.950 9049 0.6198 —0.916 6798 0.1491 
—0.938 .6728 0.7476 —0.944 9112 0.6218 —0.910 6888 0.1520 
—0.926 . 6804 0.7501 —0.936 9198 0.6246 =0-922 6979 0.1550 
=05951 6901 0.7532 —0.937 9264 0.6268 0.99% 954.7309 0.7773 
—0.940 6995 0.7563 —0.944 9355 0.6298 —0.924 7434 0.7814 
—0.933 . 7356 0.7681 —0.942 9424 0.6320 =0:937 1528 0.7845 
—0.931 . 7457 0.7714 —0.942 9709 0.6413 —0.925 7625 0.7876 
—0.941 jfool 0.7744 —0.928 . 9806 0.6445 —0.928 7716 0.7906 
+0.222 624.6542 0.0266 —0.946 9875 0.6467 —0.920 . 7809 0.7936 
—0.079 6726 0.0326 —0.940 9973 0.6499 —0.933 . 7900 0.7966 
Or 227, 6830 0.0360 —0.958 841.0046 0.6523 —0.937 1986 0.7994 
=05352 6938 0.0395 =0.953 0146 0.6556 —0.920 8073 0.8023 
—0.941 818.8972 0.4356 —0.926 860.8460 0.1294 —0.929 8177 0.8056 
—0:957 9111 0.4401 —0.909 8554 0.1326 —0.922 8264 0.8085 
—0.951 9177 0.4422 = 02929 8648 0.1356 —0.920 8365 — 0.8118 | 
—0.942 9271 0.4453 =0:912 8745 0. 1387 —0.913 8455 0.8147 
—0.944 9344 0.4477 —0.920. . 8863 0.1424 —0.909 8545 0.8177 
—0.944 9437 0.4507 —0.921 + 862.7884 0.7636 


verting the elements derived from the solution of one 
eclipse to those which should follow from the solution 
of the other eclipse involve the eccentricity itself and 
the longitude of periastron. It thus becomes important 
to know the aspect of the orbit at the time the obser- 
vations were taken. From one point of view this follows 


most conveniently from the displacement of secondary — 


eclipse from the half-period point. The previous inves- 
tigations of TX UMa by O’Keefe and Wood indicate the 
locations of secondary eclipse to be about at the ex- 
pected phases. Inspection of the shallow minimum in the 
Tucson light curves, illustrated in Fig. 2, shows that the 


TaBLE II. Magnitudes and color indices of comparison stars and TX UMa. 


Star Vv B-V Ue n 
+47°1797 +7=38-+-0002 +0"15-0™01 +0"08--0"01 10 
+47°1799 8.5640 .02 +0.124-0.01 +0.09+0.01 10 
TX UMa, 0200 8.76 +0.47 +0.01 oe 
TX UMa, 0.25 7.06 +0.01 —0.36 
TX UMa, 0.50 7.13 —0.02 —0.37 
TX UMa, 0.75 +7.06 +0.01 —0.36 
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| Taste III. Provisional yellow solution of TX UMa. (V-c) 

| any 

«+=0.6 j=82°3 -0.35 

| k=0.615 apt= 0.945 

| a,=0.270 L,= 0.196 

a,=0.166 L,= 0.804- 
-0.30 


slipse is displaced in the proper direction and that its 
ication seems not to be a function of the wavelength 
£ observation. Detailed study of the observations fails, 
owever, to convince one that the lowest point of 
clipse occurs at the predicted phase. Indeed, it seems 
s if the eclipse is displaced only about half the expected 
eparture and that the light curve varies conspicuously 
1 shape within eclipse. As a consequence, information 
erived from the analysis of secondary minimum must 
viewed with considerable reservation. More than 
kely, the eccentricity of the orbit and the longitude of 
eriastron are really the values which are expected from 
he discussions of times of primary minimum, and 
omehow secondary eclipse is distorted in shape and 
erhaps in depth. It should not be surprising, therefore, 
f inconsistency arises between the separate analyses 
or the elements from the two eclipses. For the present 
rientation of the line of apsides the effect of orbital 
ccentricity upon the eclipses creates only a small 
leparture from a circular solution. 

Each of the yellow, blue, and violet light curves was 
ndividually rectified for darkenings of 0.2, 0.4, 0.6, 
nd 0.8. Trials using the Princeton nomographs revealed 
he solution for each light curve to be similar to that 
ound by previous investigators. Attempts to refine the 
\omograph solutions were also made in each instance, 
nd for this purpose the y tables of Merrill (1950) were 
ised. By and large, these attempts were successful, and 
elected elements of the provisional yellow solution, 
eferred to the eccentric orbit, are listed in Table III. 
t should be remarked that the determination of the 
imb darkening coefficient is not strong. The fit of the 
olution suggests that a relatively large darkening gives 
he best results, and this information is substantiated by 
he coefficients predicted by Miinch and Chandrasekhar 
1949). A more detailed comparison of theory and ob- 
ervation by Hosokawa (1957a) suggests further that a 
mean coefficient of 0.6 works just as well as the indi- 
vidual ones which follow from his theory. 

Solutions for the blue and violet observations were 
ilso obtained using the y tables. In general, these 
letermine geometric and photometric elements which 
ire similar to those determined from the yellow obser- 
vations. There are, however, disturbing discordances 
which prevent exact consistency of one light solution 
with. another. As a consequence of this, it was decided 
-o test the mean solution against the theoretical Russell 
model and thus seek some clue to the sources of the 
jiscordances. 

From the depth of primary eclipse used in the pre- 
iminary solutions at each wavelength and from the 
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Fic. 2. Yellow and violet observations of TX UMa within 
secondary eclipse. The filled circles refer to observations of the 
first season and the plus signs to observations of the second season. 
The phases of the observations have been computed from the 
ephemeris of Sec. III. The theoretical location of secondary 
minimum computed from the same ephemeris is indicated by the 
arrows. 


U,B,V photometry of the comparison star it is possible 
to derive the color indices of the components of TX 
UMa. This is shown in Table IV, in which the first 
column lists the stars individually ; the second and third 
columns list the color indices of the interior, facing 
hemispheres; and the fourth and fifth columns give the 
color indices for the exterior, remote hemispheres. It 
can be seen that, in general, the interior hemispheres are 
bluer (i.e., hotter) than the outside ones. The differences 
are greater in B—V than in U—B, and they are greater 
for the cool star than for the hot star. The hot star 
seems clearly to exhibit colors associated with B stars, 
suggesting that the spectral classification of B8 is well 
established. The colors given for the cool component 
place the range of possible spectral class from F8 to G2. 
In all probability, the luminosity class of this com- 
ponent is that of a subgiant, as can be seen from the 
ratio of the radii determined from the provisional 
solution. These conclusions are in general agreement 
with those of Wood, who found the ratio of the eclipse 
depths in blue light to determine a spectral type of 
about F9 for the cool star. He suggests that the dif- 
ference between the spectroscopic classification of gF2 
and that of a late F star is significant. 


TaBLeE IV. Preliminary color indices of components of TX UMa. 


Star B=V U-—B (B—V)’ (U—B)’ 
Hot —O710 —0™40 —0708 —0™39 
Cool +0.62 +0.04 +0.70 +0.04 


TABLE V. Reflection effect in cos@ for TX UMa. 


ROBE Rae 


KOCH 


Light curve (E,/E.)R' (En/E.)x’ Air’ Ain’ Ay (En/ Ee)” (En/ Ee)” 
Yellow 0.62 0.54 —0.037 —0.042 —0.019 1.08 1.06 
Blue Lats 1.10 0.021 0.022 0.014 1.59 1.54 
Violet 1.74 185 —0.014 —0.019 —0.016 1.49 1.46 


= 


With the spectral types of the stars known more 
surely than at the beginning of the solution, it is possible 
to return to the outside-eclipse variations and compute 
the theoretical coefficient of cos@ at each wavelength. 
This coefficient is defined as 


A,=—0.40[Lir 2 (E./En)— Len? (En/E-) | sini. 


In this equation r refers to the mean radius of a com- 
ponent; L refers to the luminosity of its interior hemi- 
sphere; and E to the luminous efficiency of its surface 
as it reradiates the light that falls upon it from the 
companion star. The L’s and 7’s follow from the 
provisional solutions, the yellow values being given in 
Table III, and the luminous efficiencies are given by 
theory. At present, two theories for the computation _of 
E are available. There is, first of all, that developed by 
Russell and Pismis and quoted in Russell and Merrill 
(1952c) which depends upon graybody radiation for the 
star. There is, in addition, the theory developed by 
Hosokawa (1957b), which takes account of the de- 
partures of the stellar radiation from graybody char- 
acter. The computed ratios of the luminous efficiencies 
are shown in Table V, where the subscripts H and R 
refer to the values determined by Hosokawa’s and 
Russell’s theories, respectively. From these computed 
ratios the theoretical values of A; have been calculated. 
They are shown in columns four and five of Table V 
and are contrasted with the observed values of A; 
rounded off to three decimals. It can be seen that the 
theoretical values of the differential reflection coefficient 
differ significantly from the observed values and that 
the discordance becomes worse as one goes tq longer 
wavelengths. This is a situation which is, at least in 
part, familiar to many observers of eclipsing stars. For 
instance, a similar situation has been reported for Algol 
by Hosokawa (1960). It is also well known that in par- 
ticularly close binaries, such as those of the W UMa 
type, the coefficient of reflection is frequently positive 
in contradiction to the theoretically negative value. The 
situation in TX Uma, at least for yellow and blue light, 
is similar in that the observed values are more positive 
than the theoretical ones. The main assumptions which 
underlie the theory of reflection as it is presently applied 
are that the stars are spheres and that each is in perfect 
radiative equilibrium. There are, however, additional 
secondary effects which can be treated only} In approxi- 
mation. Among these are the quantitative allowances 
for the finite separation and sizes of the components 
together with penumbral effects of illumination. No 
improvement in the treatment of these being expected 


at present, it seems that the observed departures should 
be attributed to a grossly incorrect solution or tq 
violations of at least one of the two main assumptions} 

Fortunately, the terms which make up the expressior} 
for A, are at least partly separable so that they can be 
checked by independent means. If we suppose first of 
all that it is the relative luminosities of the stars which 
are in error, it is possible to juggle these so as to 
reproduce the observed reflection coefficients. This 
occurs in blue light when L;, is approximately 0.7 and 
L, is 0.3. This explanation, however, cannot be the 
correct one, for outside eclipse we see the absorption 
lines of only the hot star. It is possible that the radii of 
the stars given by the preliminary solution are in error. 
These, however, agree satisfactorily with previous 
determinations. 

It is difficult to see how a departure from spherical 
figure can decrease the amount of light which is reflected 
toward the observer. As the stars are distorted by their 
rotations and mutual gravitational attraction, the 
interior hemisphere of each is elongated more than the 
exterior side. This interior hemisphere, being the main 
seat of reflected radiation, will therefore be larger than 
in an undistorted sphere, and the total amount of 
reflected light, if it obeys the distribution law given by 
Milne (1926) and others, should be greater than that 
reflected from a spherical star. More than likely, there 
are appreciable departures from radiative equilibrium 
and theoretical luminous efficiency within the atmos- 
pheres of the stars. In particular, the cool star of larger 
size and lower surface gravity is suspect. A different 
conclusion has been expressed by Hosokawa (1957b), 
who suggested that the source of this same difficulty, 
encountered in his statistical study, can be traced to 
incorrect values of E for B stars. This conclusion, how- 
ever, is not consistent with the situation in TX UMa, 
for the discordance in A; becomes greater as the lumi- 
nosity of the cool star contributes a larger fraction of 
the light of the binary. It seems more likely, therefore, 
that the efficiencies which are incorrect are those asso- 
ciated with the cool component. This being so, it is 
possible to determine the ratio, E,/E., from the 
observed luminosities and radii of the stars given by the 
preliminary solutions and the observed values of the 
differential reflection coefficient. These forced ratios of 
luminous efficiency are given in the last column of 
Table V. If we take at face value the theoretical effi- 
ciences tabulated for a B8 star, then those determined 
for the cool component by the observed reflection coef- 
ficients are ones associated with a G or K star. To sum 


p, the view presented asserts that a preliminary 
~ lution based upon either graybody or nongraybody 
 neory will give an informative comparison of the 
neoretical and observed reflection coefficients. Should 
‘nese be found to differ, the ratio of the luminous 
| fliciencies as a function of wavelength may be deter- 
uined from the observed reflection effect. It follows 
hat the luminous efficiency is probably a function of 
he sizes and the separation of the stars and should be 
jetermined individually for each eclipsing binary. This 
ppears to be not too difficult to accept, for the theo- 
etical values computed according to either one of the 
vheories pertain to single, isolated stars. 

| It is now necessary to compare the observed and 
\ heoretically expected values of the coefficient of cos26, 
/onventionally A», in the observed light curves. This 
oefficient is predominantly determined by the ellip- 
‘icities of figure of the stars, but there is, in addition, a 
small term contributed by the reflection effect. The 
‘omplete expression is given, for instance, by Hosokawa 
‘(1957a) in his Eq. (17). For the purposes of the present 
“nvestigation, his final term standing for the small 
ourth-order distortion of figure is neglected. The mean 
‘photometric ellipticities follow from the observed 
outside-eclipse variations, and these have been used for 
the preliminary solutions of the eclipses themselves. As 
has been pointed out by Dr. J. E. Merrill in private 
conversation, it is almost certain that the two com- 
ponents of TX UMa differ significantly in their figures. 
This assertion can be examined by the tests given by 
Russell and Merrill (1952d). At the University of 
Pennsylvania, Merrill and Reuning have calculated, 
but not yet published, the radii in four directions for 
the m=3 polytropes of Chandrasekhar (1933). From 
the mass ratio given by Pearce (1933) and Wood and 
the radii given by Merrill and Reuning it is possible by 
interpolation and iteration to recover the radius of the 
undistorted polytrope and then compute the two 
equatorial radii in the direction of the line of centers, 
the equatorial radius in the plane perpendicular to the 
line of centers, and the polar radius. It is found that for 
each component of TX UMa the radii in the direction 
of the line of centers are very nearly equal. Finally, it 
is possible to compute the equatorial and polar ellip- 
ticities of figure, customarily represented as z. For the 
hot star the two ellipticities are each equal to 0.011. For 
the cool star the equatorial ellipticity is 0.187, while the 
polar ellipticity is 0.093. It is thus confirmed that the 
two stars are distinctly different in figure, and that the 
cool star is distorted from an ellipsoid of rotation. At 
this stage of the computation it seemed permissable to 
work with a mean value of the ellipticity of the cool 
star, and this was taken as 0.139. : 

In order to compute the expected value of the coef- 
ficient of cos26, it is necessary to make some assumption 
about the ratio of photometric ellipticity to geometric 
ellipticity. This is a quantity, little known today, which 
depends upon the gravity and limb darkenings of a star. 


TX UMa 
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TABLE VI. Reflection and ellipticity effect in 
cos26 for TX UMa. 

Light curve A, Ag!” A» Alt? 
Yellow —0.007 —0.014 —0.014 —0.014 
Blue 0.007 0.010 0.006 0.014 
Violet —0.006 —0.006 —0.006 —0.010 


Investigation by Hosokawa (1957a) again shows that, 
in general, one can do little better than choose a value 
for this ratio appropriate to a limb darkening of 0.6. 
From these quantities and the ratios of the luminous 
efficiencies, the values of A» have been computed and 
are listed in Table VI as a function of wavelength. For 
this computation the disparities between the luminous 
efficiencies of Russell and Hosokawa are of no sig- 
nificance because the ratios are multiplied by a small 
number. These theoretical values are given as Ao’ in 
Table VI, while the values computed with the forced 
luminous efficiencies are given as A’. The observed 
values are listed as A>. It can be seen that the change in 
the luminous efficiencies for the violet curve produced 
no change in the computed ellipticity effect. For the 
blue curve the agreement with the new luminous 
efficiencies is slightly worse than it was previously, but 
for the yellow curve the agreement is now strikingly 
better. This is taken as evidence that forcing the ratios 
of luminous efficiency in order to fit the individual 
binary is a change in the proper direction. 

The elements of the system being provisionally 
known, it is possible to compute the theoretically 
expected coefficients of the higher-order terms in the 
reflection and ellipticity effects. These depend upon the 
luminosities of the individual stars and their radii and 
the mass ratio. The coefficient of the third-order term 
for TX UMa is computed to be about —0.001 though 
there is a slight variation with wavelength as a con- 
sequence of the dependence of the luminosities on wave- 
length. The fourth harmonic of the ellipticity and reflec- 
tion effects is computed to be 0.000. 

The attitude having been taken that the ratios of the 
luminous efficiencies of the components of TX UMa 
have been improved, it was decided to re-rectify the 
individual observations, and accordingly the outside- 
eclipse variations were inspected once again. It was 
concluded that any improvements to be made in the 
coefficients were small and would not affect the solutions 
significantly. At this time it was found useful to change 
by a small amount the normalizing magnitude dif- 
ferences to be added to the blue and violet observations 
and thus change these outside-eclipse variations by a 
factor close to unity. These new variations are given as 


£(b) =0.9852—0.0139 cosd—0.0056 cos28, 
£(u) =0.9846—0.0163 cosd—0.0059 cos26. 


The new coefficients for rectification are given in 
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Taste VII. Improved rectification coefficients for TX UMa. 


Light curve Co Cy Cy Nz 
Yellow +0.0185 +0.0192 +0.0059 +0.0812 
Blue 0.0135 0.0139 0.0043 0.0407 
Violet +0.0143 +0.0163 +0.0047 +0.0432 


Table VII. As anticipated, changes in the rectified ob- 
servations are small. 

A determinate, consistent solution from the com- 
bination of both eclipses was now sought. It was found 
that, while at each wavelength there is indeed a deter- 
minate solution defined by the depth of both eclipses 
and the shape of primary, the solutions so discovered 
are not identical and cannot be made identical within 
the observational error by changing the darkening coef- 
ficient. At least some of these difficulties can be traced 
to the shape of primary eclipse, which varies with 
wavelength. For a second trial it was decided to seek the 
best common solution defined only by the depth lines 
for all wavelengths. Again, it was found that no con- 
sistent, recognizable, partial solution results so that the 
observed depth of at least one of the eclipses is sensibly 
incorrect in at least one wavelength. Since the depths 
of primary eclipse give B—V and U—B color indices 
that are reasonably appropriate for a B8 star, it seems 
likely that a large percentage of the discrepancies is 
to be traced to incorrect depths of secondary. In order 
to test this, the color indices for the cool component 
were derived from the depths of secondary eclipse. 
These were found to be 072 to 03 bluer than should 
be the case for a late F-type star. As a result, it seems 
certain that the observed depth of secondary eclipse is 
significantly incorrect in at least two wavelengths and 
one would like to know which of the observed depths is 
most consistent with primary eclipse. This it is also 
possible to find, and three trials determined that the 
depth of secondary eclipse observed in violet light is 
most nearly correct. The departures of the secondary 
eclipse depths from what they should be seem to be 
correlated with the departures of the reflection coef- 
ficients from the theoretical values. This occurs in the 
sense that the observed reflection coefficients make the 
depths of secondary too shallow in yellow and blue light. 

So many difficulties having been detailed already, it 
may be useful to summarize them at this point. First, 
though the depths of primary eclipse appear reasonably 
accordant, it seems certain that there are small dis- 
parities in the shapes which make it virtually impossible 
to determine a unique, consistent solution. Second, 
between eclipses the light variations differ from those 


computed on the basis of the Russell model. Agreement 


may be somewhat improved if empirically determined 
luminous efficiencies are used instead of those which 
follow from blackbody, graybody, and nongraybody 


ROBERT H. KOCH 


theories. Third, the depth of secondary eclipse ij) 
certainly not what theory predicts in yellow and bla 
light, but agreement is satisfactory in violet light}, 
Finally, the shapes of the branches of secondary eclips}j) 
suggest strongly that light variations may be promineny, 
at this phase in the light curve and that they too may) 
be a function of wavelength. These discrepancies beinj\) 
understood to inhibit the determination of any solution ) 
it seems the most logical procedure to rely on the depths} 
and shapes of primary eclipse as far as they can bi 
brought to consistency and on the observed depth oiff 
the violet secondary eclipse. 


characteristics traceable to an atmospheric eclipse. His}) 


conclusion, while it cannot be validated for the time of 


absorption capable of distorting the shapes and the 


depths of at least one of the eclipses. It was finally'# 


decided to treat the yellow observations in primary for 
a body eclipse involving sharp-edged stars but to} 
permit a sharp-edged, extensive atmosphere to distort | 
this model for the blue and violet observations. The | 
atmosphere is conceived to absorb selectively, 
may also absorb and scatter in a nonselective fashion. 
Though this cannot be proved to be the correct attack, 


it is certainly true that the absorption lines normally | 
present in stellar spectra are more numerous toward | 
shorter wavelengths, so that atmospheric effects should © 
be enhanced in passing from yellow to blue to violet | 
observations. To make such a solution it is necessary | 


to depart from the Russell model in favor of one similar 
to that expounded by Kron and Gordon (1950) in 
their study of V444 Cyg. In their model both stars are 
sharp-edged, but one is surrounded by at least a less 
dense shell. In projection the rectified model exhibits 


A 
but it | 


: 


| 


7 


two circular stellar disks which may be darkened to the ] 
limb in the ordinary manner, and the component pos- — 
sessing the extensive atmosphere is surrounded by an — 


annulus which is concentric with the star and touches 
it at its photosphere. The equations which permit the 
solution of primary eclipse, in the case of TX UMa the 
partially atmospheric. occultation eclipse, are given 
below for equal, unspecified darkenings in yellow and 
blue. In these equations the unprimed symbols refer to 
characteristics of the photospheres of the components. 
The singly primed symbols refer to the disk described 
by Kron and Gordon, and the doubly primed symbols 
refer to the projected atmospheric annulus. The equa- 


tions for yellow and violet are similar, and the extension 


to dissimilar darkenings is performed by the usual 
technique. The arrows indicate the sequence of calcu- 
lating the succeeding quantities. 


‘ j 

bli (a,k,p)°°= (a,k,p) 0° 70, 

i (5=0', k’)s— pr ay’ %, kh’ =k", 
IB ay’ 0c = q,’ 08 — Gy, °°, 


| (1—£),°+ (1-0”),°°=d, 
‘| (1—€0)0°°= (a0,6°°) (Las), 
 A-f")5"= Af) aoa” *) (Las) fe 


[ere d equals the observed depth of primary eclipse in 
lue light at any phase. 
| The solution is begun by assuming provisional 
spouts of atmospheric absorption in blue and violet 
‘nd is continued by successive trial values for the 
‘mgitude of external tangency, (1—()’’)*, and a,’’. It 
yas decided that physical reality is best served by 
i, ssuming the limb darkening coefficient in blue and 
ec to be not less than that chosen in yellow. The 
diest yellow solution is tabulated in Table VIII. For 
/.6 limb darkening it proved impossible to fit the blue 
nd violet observations with a consistent body and 
».tmospheric eclipse as well as the observations sug- 
tested might be done. Trial solutions for both blue and 
riolet light were next computed assuming 0.6 limb 
)larkening for the photospheres and 0.8 darkening for 
the atmospheres. These proved as inferior as the 
yrevious 0.6 solutions. Solutions were obtained, how- 
‘ver, when the blue and violet observations were 
‘reated for 0.8 darkening both for the disks and for the 
itmospheres; refined through several trials, these also 
appear in Table VIII. In those characteristics which 
tefer to the sharp-edged photospheres the solutions at 
all wavelengths are now consistent. The light curves 
were analyzed on a scale 1 cm to 0.01 in rectified light. 
Reduced to the size of a printed figure the differences 
between unsatisfactory and satisfactory solutions 
cannot be distinguished. For a single observation the 
probable residual from the final solution is +0™010 
in yellow, +0™009 in blue, and +0008 in violet light. 


fi 


Tare VIII. Final solution for primary eclipse of TX UMa. 


k= 0.615 V,/=+7732 
: a,= 0.270 V./=+9.00 
b,= 0.267 (B—V),/=—0.09 
as= 0.166 (B—V),’=+0.58 
b,= 0.164 (U—B),’=—0.41 
j=82°2 (U—B),’=+0.21 
y b u 
(1—£:)% 0.7570 0.8296 0.8647 
x 0.6 0.8 0.8 
a 0.945 0.953 0.953 
Lr=Ls 0.801 0.871 0.908 
L.=L, 0.199 0.129 0.092 
3 0.031 0.013 0.014 
Ly’ 0.796 0.866 0.905 
Lr! 0.156 0.098 0.057 
Fiat Je das 13:5 23.4 41.8 
(1—£9"") PF oes 0.0014 0.0033 
eae th 0.8 0.8 
ay’ °° 0.031 0.031 
7, 0.948 0.883 
"1 0.285 0.285 


| Tx! UMa 
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The large residual in yellow is caused, in large part, by 
the scatter of the observations at the bottom of eclipse. 

The final values of (£;/E,)/” and A,’” are given in 
Tables V and VI, respectively. The former have 
changed little. The latter are now too negative in blue 
and violet, a situation that can be traced to the necessity 
for 0.8 limb darkening at these wavelengths. 


VII. DISCUSSION 


It is highly desirable to comment upon the deter- 
minancy of the elements listed in Table VIII. While 
extensive treatment of the yellow observations demon- 
strated that darkenings of 0.2 and 0.8 were not the 
best that could be desired, it proved impossible to dis- 
criminate between values of 0.4 and 0.6. The solution 
thus rests on the assumption that the true coefficient is 
0.6. The difference of darkening between the yellow and 
blue or violet sets of observations is undeniable, but 
it has not been proved that this difference would exist 
if both components were normal, well-separated stars. 
In other words, it may be that the observed difference 
is an artifact of the complicating, atmospheric effects 
that we see in primary eclipse. As an example of the 
determinancy of the yellow solution of primary itself, 
it can be remarked that the following solutions on the 
assumption of a 0.6 limb darkening were found to be 
distinctly inferior: k=0.602, ao°=0.955; k=0.620, 
a1°°=0.900; and k=0.635, ac°°=0.895. These were not 
the only solutions so tested, but they are the ones closest 
to the final solution. There is a certain amount of latitude 
in determining the best atmospheric solutions in the 
blue and violet. Trial solutions serve to indicate, how- 
ever, that a,’’, for instance, is not 0.277 or 0.293. It is 
legitimate to conclude then that this particular element 
is known to an accuracy of about 3%. There is no 
physical constraint which necessitates that a,’’ be the 
same in blue and violet light, and probably a difference 
could be detected with more accurate observations. The 
eclipse observations of the second season, while few in 
number and covering but the top portion of one of the 
shoulders of the light curve, demonstrate that primary 
minimum varies in shape. Furthermore, the departures 
of the observations of this season from those of the 
first one are greater in blue and violet than they are in 
yellow. One is tempted to suspect, therefore, continuing 
activity associated with the cool component which can 
be treated to a reasonable degree of accuracy as growth 
and decay of an extensive envelope. It is quite likely 
that there is a certain amount of atmospheric absorption 
present in the yellow light curve. In the absence of 
extensive spectroscopic observations in the visual region 
this feature remains unknown. 

From the limiting solutions given by Wood the 
radius of the cool star seems reasonably well known, so 
that the depth of the atmosphere can be calculated. 
This is 1.53 10° km. The line of sight from the observer 
to the disk of the hot star traverses at most 1.86 10° 
km through the extensive atmosphere of the cool star. 
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It was thought at first that these data, together with 
plausible values of the mean molecular weight and 
temperature in the atmosphere, could provide some 
information about the density gradient with atmos- 
pheric height. Because of the small number of the 
observations and the inevitable inaccuracies of these, 
this expectation has not been borne out. 

It should finally be remarked that the observations 
were tested against the atmospheric tables of Linnell 
(Astrophys. J. to be published). These tables permit 
the application of five different density gradients in an 
extensive atmosphere. It was found that none of these 
could be made to represent the observations of TX 
UMa, thus indicating the presence of a photosphere for 
the cool component. 


VIII. SECONDARY ECLIPSE 


Since the use of the observed coefficients of reflection 
gave disparate colors for the cool star, it seemed useful 
to rectify at least secondary eclipse according to the 
average theoretical values of A; given in Table V. This 
treatment of these observations by the admission of a 
complication term in cos# is not completely successful 
for two reasons. First, the coefficient of cos2@ to be used 
in the rectification for ellipticity in blue light comes 
out to be +0.0005 and thus of the wrong algebraic 
sign. Second, the shape of each branch of the minimum 
is still inconsistent with the solution from primary 
eclipse. There is, however, a strong indication that this 
piecemeal rectification is the correct procedure to 
follow, for the colors of the cool star as they can now 
be derived from the rectified depths of secondary are 
those of a late F-type star as were the colors derived 
from primary eclipse. 

Though this ad hoc method of treating secondary 
eclipse can do nothing about the strongly asymmetric 
minimum at each wavelength, it does suggest a reason- 
ably coherent model for TX UMa. A partially absorb- 
ing, re-emitting, and fluorescing atmosphere surrounds 
the cool star. Its presence and the activity within it 
create departures from the theoretical reflection effect 
and the ellipticity effect insofar as the latter depends on 
the observed reflection and limb darkening. At second- 
ary eclipse the uneclipsed fraction of the atmosphere is 
greater than the eclipsed portion, causing (1—¢)* to be 
too small. The variable extent and density of the 
atmosphere cause changes in the depths and shapes of 
the eclipses, but the changes | in Soy eclipse are 
proportionately more pronounced. 


Ix. SUMMARY 


It should be fairly clear that the solution given in 
Table VIII cannot be considered the definitive one for 


fie OCH 


TX UMa, but it does emphasize a few points i | 
bear on the treatment of all eclipsing variables. Par} 
ticularly for those in which one component nearly fill 
the inner Lagrangian surface, observation at only oni} 
wavelength offers no assurance that correct element a 
can be determined. This forces us to view with con} 
siderable reservation orbital solutions which rest ot} 
even well-covered visual or photographic light curves} 
At best, these provide an instantaneous photometri¢ 
cation: | 
Finally, if the combination of body and atmospherid 
eclipses is shown to have extensive application, it might 
be profitable to formalize the Kron model more ele: 
gantly and compute tables for it. | 
/ 
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The internal constitutions of Uranus and Neptune have been investigated on the basis of the assumption 
that solid ammonium is a major constituent of these planets. An equation of state found by Bernal and 
Massey for this substance has been employed together with the equation of hydrostatic equilibrium to find 
the internal density and mass distributions. The presence of undetermined, more dense materials has been 
allowed for in a phenomenological manner and it has been found that the ammonium contents are 84 and 
747% for Uranus and Neptune, respectively. Both planets appear to have dense cores of heavy material of 


I. INTRODUCTION 


J T is the purpose of this paper to report the results 
4% of an investigation of the internal structure of 

ranus and Neptune based upon the hypothesis that 
metallic ammonium is a major constituent of these 
. This possibility has been pointed out by 


showed that Uranus and Neptune cannot be composed 
primarily of solid hydrogen, as are Jupiter and Saturn, 
ut are likely to contain large amounts of water, 
ammonia, or methane. Bernal and Massey (1954) have 
shown that ammonia (VH;) and hydrogen will combine 
to form stable ammonium metal (NH,) at pressures 
less than those which prevail throughout most of the 
interiors of Uranus and Neptune, and have thus given 
further support to Ramsey’s suggestion. The former 
workers have also calculated by quantum-mechanical 
methods the 0°K isotherm of the metal and a quanti- 
tative study of these planets may be made. 

The only current models of the two smaller major 
planets are those of Wildt (1947), who has assumed that 
they were composed of three layers of material, each 
layer or shell having a uniform density. When reasona- 
ble values were assigned to these densities, the radii of 
discontinuity could be calculated by requiring the 


TABLE I. Pressure-density relation for metallic 
ammonium at 0°K. 


Pressure Density 

(102 dyne/cm?) (g/cm) 
0 00 0.85 
0.04 28 
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*Part of this work was performed while the author was a 
summer visitor with the Theoretical Division, Goddard Space 
Flight Center, National Aeronautics and Space Administration, 
Silver Spring, Maryland. 


models to reproduce the correct total planetary mass 
and moment of inertia. The latter quantity may be 
found by the Radau-Darwin formula from the observed 
ellipticity and rotational velocity of the planets. That 
such models may lead to erroneous conclusions regard- 
ing the matter of which the planets are composed was 
shown by Ramsey (1951, 1952) and by De Marcus 
(1958) for the cases of Jupiter and Saturn. It therefore 
seemed desirable to study the properties of more realis- 
tic models of Neptune and Uranus. 


II, PROPERTIES OF METALLIC AMMONIUM 
AND EMPIRICAL DATA 


In their calculation, Bernal and Massey (1954) first 
obtained a self-consistent field for the ammonium mole- 
cule by averaging the tetrahedral distribution of 
protons over a sphere of radius 79. The appropriate 
value of 7) was obtained from a calculation of the total 
energy of the tetrahedral configuration as a function of 
this radius. The energy of the ammonium metal was 
then calculated by dividing the metal into spheres of 
radius a centered on each nitrogen nucleus. By placing 
the proper requirements on the wave functions of the 
metallic electron within and on the boundaries of each 
sphere, the energy per nitrogen atom E£ was found as a 
function of a by an iterative procedure. The equation 
of state of the metal at 0°K was found from the relation 


p= (3/4r)dE/d(a'). 


Their results are given in Table I. These results are 
applicable to studies of the planets since it is unlikely 
that their internal temperatures attain values which are 
high enough to modify appreciably the 0°K isotherm.. 

In Table II are given the results for the equation of 
state of solid hydrogen found by De Marcus (1958) and 
which will be used later in the present work. 

The relevant empirical data for Uranus and Neptune 
are given in Table III. These data are taken from 
De Marcus (1959), except for the value of J/M Ro? for 
Uranus. The value of this important quantity is not 
well known and has been taken to be 0.24 as was done 
by Wildt (1947). The true value is probably close to 
this number. 
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Taste II. Pressure-density relation for solid hydrogen at 0°K. 


Pressure Density 
(102 dyne/cm?) (g/cm*) 
0.00 0.089 
0.40 0.399 
0.60 0.455 
0.80 0.780 
1.00 0.82 
2.00 1.01 
3.00 1.16 
4.00 1.29 
5.00 1.40 
10.00 1.87 


Il. CALCULATION 


Under the assumption of hydrostatic equilibrium in 
the interior of a cold, rotating body, any pressure p 
defines an oblate isobaric surface which is approxi- 
mately a spheroid. If 7 is the radius of a sphere which 
has the same volume as the spheroid and if p is the 
density on the surface of this sphere, then, for small 
rotational velocities w the equation 


dp/dr=—p(r)Lg(r)—3o"r ] (1) 


must hold. Here, g(r) is the local acceleration due to 
gravity and is given by 


g(r) = (4nG/r") A p(x)ade, (2) 


where G is the gravitational constant. It is understood 
that (1) is to be solved together with the equation of 
state connecting the pressure and density of the matter 
of which the planet is hypothesized to be composed. 
The result of such a solution is the knowledge of the 
density and pressure as functions of r. The mass dis- 
tribution may be found from the relation 


M(r)=(7°/G)g(r), 


while the moment of inertia may be calculated from 


(3a) 


T= (8n/3) i Rete (3b) 


It has been found to be convenient to perform the 
numerical integration of (1) by beginning at the outer 
planetary radius and carrying the integration inward 
to the center. The appropriate starting conditions were 


Taste III. Empirical data of Uranus and Neptune. 


Mean Mean 
Mass radius density 
Planet =M)(10?7 g) =Ro(km) (g/cm’) T/M Re 
Uranus 86.9 23 700 1.56 0.24 
Neptune 102.8 21 500 2.47 0.30 
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at r= Ro. (4h 


p=p (p=0) 


r=0.05 Ro and the remaining mass was ascribed to ¢ 


4 


heavy central “particle” of constant density. Thi 
occurrence of such a central “particle,” or core, may bq 
explained qualitatively as due to the settling undej 
gravity of the heavy elements in the planet since thd 
time of its formation. | 
Since it was desired to allow for the presence 0 
materials of higher density than ammonium, the densit} 
| 


function used in (1) was taken to be | 
e(r)=pi(r) LI t+7(1—7/Ro) J, (5) 


where py; is the density of the ammonium and ¥ is con- 
stant. This form is also in accord with the idea of gravi- 


assuming presence of ammonium plus heavier materials (y=0.0 


TaBLE IV. Calculated properties of the interior of Uranus 
I/M RP? =0.32). j 


p 
s=r/Ro (g/cm?) M/Mo 

1.00 0.84 1.000 0.000 | 
0.90 1.07 0.834 , 0.227 ‘ 
0.80 1.43 0.661 0.531 : 
0.70 1.74 0.489 0.915 | 
0.60 2.01 0.336 135 
0.50 221 0.212 1.80 
0.40 2.36 0.123 2.24 
0.30 2.48 0.065 2.67 
0.20 ZO Sica 0.034 Sai) 
0.10 2.70 0.023 3.97 
0.05 Des ff) 0.021 5.40 


tational settling of the heavier elements. The constant 
y has been treated as an adjustable parameter which 
was chosen so as to give the best agreement between 
the computed and empirically determined values of I 
and to imply the smallest core mass. 

The important characteristics of the best models ob- 
tained in this way for Uranus and Neptune are given 
in Tables IV and V. The results for Neptune are thought 
to be acceptable but the calculated value of J is too 
large by 30% for Uranus. ; 

In order to obtain an acceptable model of Uranus, it 
was necessary to assume the presence of some material 
of smaller density than ammonium, at least in the outer 
regions of the planet. Since its equation of state is well 
known, this material was taken to be solid hydrogen. 
If ideal mixing is assumed, the density in the presence 
of both ammonium and hydrogen is given by 


1/p= (1—w)/prt+-w/p», (6) 


where pz: is the density and w is the mass fraction of solid 
hydrogen. The parameter w has been treated as an ad- 


| 


: 
jstable constant and it has been further assumed that 
> radii less than 71, w became zero. Several values of 
@iand r; permitted acceptable models of Uranus, the 
Pst satisfactory of which is given in Table VI. 

| Poor results have been obtained by the application of 
tese assumptions to Neptune. 


| IV. CONCLUSIONS 


i It is concluded that acceptable models of Uranus and 
‘eptune may be constructed in which metallic ammo- 
um is a major constituent. The best case for Neptune 
i rresponds to an ammonium content of 74% by mass, 
e remaining material being composed of heavier ele- 
pats That for Uranus corresponds to an ammonium 
patent of 84%, a hydrogen content of 2%, and a 
~avier element content of 14%. It is not considered 
‘gnificant that Uranus appears to have retained some 
eiozen while Neptune, the more massive body, has 
a It is thought that allowance for the simultaneous 


-esence of materials both more and less dense than 
“nmonium within the interior of Neptune would lead 
i an acceptable model of that planet in which some 
ydrogen would occur. However, due to the present 
ie of information concerning the pressure-density 
‘lationships of the heavier elements at high pressures, 
| would seem that such a model would involve an un- 
ortunate number of adjustable parameters. 
iit may be noted that the masses of the cores of the 
{ 


| Taste V. Calculated properties of the interior of Neptune 
Tio presence of ammonium plus heavier materials (y=1.00, 
iM oRo?=0.30). 


| p p 
| s=r/Ro (g/cm?) M/Mo (10 dyne/cm?) 
. 
Hf 1.00 0.85 1.000 0.000 
i 0.90 1.34 0.883 0.352 
— 0.80 2.14 0.732 0.970 
0.70 292 0.558 1.89 
0.60 3.58 0.391 3.04 
| 0.50 4.08 0.249 4.31 
0.40 4.35 0.144 5.58 
0.30 4.62 “0.077 6.81 
0.20 4.90 0.040 8.18 
0.10 5:17 0.026 10.6 
0.05 5.30 0.024 14.8 
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Taste VI. Calculated properties of the interior of Uranus 
assuming presence of some solid hydrogen (w=4, 71=0.90Ro, 
I/M Ro? =0.26). 


p p 
s=7r/Ro (g/cm?) M/M, (102 dyne/cm?) 

1.00 0.22 1.000 0.000 
0.90 0.94 0.935 0.097 
0.80 1.29 0.783 0.407 
0.70 1.70 0.620 0.853 
0.60 2.05 0.467 1.43 
0.50 Di oD 0.339 2.13 
0.40 2.54 0.244 2.97 
0.30 2.70 0.181 4.07 
0.20 DetZ 0.148 5.84 
0.10 2.76 0.136 10.5 
0.05 2.82 0.134 19.5 


best models given imply that, in both cases, the central 
concentration of heavy elements for these planets is 
approximately twice the mass of the earth. 

Methane is known to be present in the atmosphere of 
the major planets and it would therefore seem desirable 
to include it in any study of the interiors of the same 
planets. For this purpose, a calculation of its equation 
of state is currently under way. 
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The orbit of periodic comet Oterma has been integrated on an IBM 650 over the forthcoming very close 
approach to Jupiter. An ephemeris for 1962-66 is provided and so are osculating values of the orbital 


elements at intervals during this period. 


ERIODIC comet Oterma was discovered at the 

Turku Observatory on 1943 March 27, and pre- 
discovery images were later found on plates taken there 
in February 1942. Since then, on account of the low 
eccentricity of its orbit, the comet has been under 
observation at opposition each year. 

Dr. Paul Herget, using the IBM 650 electronic com- 
puter at the University of Cincinnati, has applied a 
differential correction to his earlier orbit (Herget 1947) 
incorporating 222 observations made mainly at the 
Turku, Lick, Yerkes, and McDonald observatories 
between 1942 February 17 and 1959 June 3. Perturba- 
tions by the earth, Jupiter, and Saturn were included. 
The resulting orbital elements, together with osculating 
elements for the 1950 and 1958 perihelion passages, are 
given in Table I. 

That the orbit is to be drastically changed during a 
prolonged close approach to Jupiter between 1961 and 
1965 has already been noted (Herget 1947; Oterma 
1958), and it was also obvious that the programs which 
are normally used at the Minor Planet Center in Cin- 
cinnati to compute perturbations by a variation-of- 
elements method would therefore be unsuitable in this 
case. Hence Cowell’s method has been programmed and 
the integration extended over the close approach, using 
the IBM 650 at Yale University. Perturbations by all 
the planets, Venus to Saturn, were included. An ephem- 
eris for the period 1962-66 is given in Table II. For 
1961 an ephemeris has been provided (Oterma 1960), 
and it is unlikely that the comet will be observable 
after early 1966 until (perhaps) shortly before the next 
perihelion-passage in 1983. The magnitudes in the 
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accompanying ephemeris have been calculated from | 
mag.= 9.5+10 logr-+5 logp, 


| 

i} 
but the comet may well be a magnitude fainter thal 
this. 

Table III gives osculating values of the orbit 
elements and also the distance (o,) of the comet fror! 
Jupiter for various dates. The results agree very closel} 
with those given by Miss Oterma (1958). It is ¢ 
interest to note that the node regresses almost 100° i 
one interval of 40 days and that for a period of sever. 
months it actually advances, that there is no addition: 
passage through perihelion, and that the size of th 
orbit is eee increased. Tisserand’s criterion (a) ha 
been evaluated for two dates, one before and one afte: 
the period of the most violent perturbations by Jupit 
It was considered in the form 


é) }} 


1 Ae 
=e La cosl, 


a a’? 

where the inclination J is referred to Jupiter’s orbit, anc 
where a’ is the radius of Jupiter’s orbit (assumec 
circular). It was found that a=0.583 in March 196 
and a=0.576 in October 1966. The difference is no 
large in view of the approximations involved, and really 
rather small when one considers that the respectivs 
values of the first term alone are 0.251 and 0.139. 

It is unfortunate that the comet is in conjunction witl 
the sun around the time of its closest approach te 
Jupiter, and that afterwards it will be so extremely 
faint; but it is to be hoped that observations can stil 


TABLE I. Orbital elements. 


1943 Oct. 3.0 
1942 Aug. 21.7241 
iss 47} 


Epoch (ET) 
T (ET) 
Bhi 0.0 


Q 

Z 

q 

€ 

a 3.961685 
n° 0.1249925 
P (years) 7.89 


1950 July 18.0 
1950 July 16.4856 


1958 June 6.0 
1958 Lee He 5045 


35427913 35428 

155.1305 1552 i100 
3.9886 3.9921 
3.404882 3.387831 
0.142763 0.144497 
3.971926 3.960045 
0.1245095 0.1250702 
7.92 7.88 
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TABLE II. Ephemeris for 1962-66. 


a (1950.0) 6 Ps r Mag. | 0“ ET a (1950.0) 6 A r Mag. 
1964 
22803™55  —11°26'8 5.504 4.639 19.9 | Feb. 5  Om59™60 +4 7°11/4 
22 12.50 10 41.1 15 1 06.26 TiS Ae Os TSM 224) 62075 
22 21.00 956.8 5.343 4.650 19.8 25 113,52 8 40.5 
22 28.96 9 14.7 Mar. 6 LUM e A EDM) Ey Cv Sno eyen ls 
22 36.30 B56) 5,120. 4.661. 19.7 
22 42.93 8 00.3 : : 1964-65 
22 48. 29. 4.852 4.671 19.6 | June 14 DSM DT Ay 2ieie 60/40 513108 207 
22 53.65 7 04.7 24 256.05 18 04.6 
22 57.54 646.0 4.559 4.682 19.5 | July 4 303.51 1839.1 5.874 5.326 20.6 
23 00.31 6 34.6 14 3 10.38 19 10.3 
23 01.87 631.0 4.267 4.692 19.4 24 3 16.56 19 38.0 52623 5.3449 20.5 
23 02.15 6 35.8 Aug. 3 3 21.95 20 02.2 
23 01.12 649.2 4.008 4.703 19.2 13 3 26.41 D297, — 5.344 52361 204 
22 58.80 7 10.8 23 3 29.83 20 39.3 
22 55.29 109-8 253.817) 4.714 19.1 | gent. 2 3:32:12 205100 95.053) 5.380 2028 
22 50.79 8 14.8 12 3 33.19 21 00.2 
22 45.58 853.7 3.725 4.725 19.1 22 3 33.00 S104. 1 4.799 ~ 5,308. - 20:2 
22 40.02 9 33.8 Oct. 2 3 31.55 21 03.5 
22 29.49 10 47.7 22 3 25.32 20 49.0 
22 25.27 1116.8 3.896 4.749 19.2 | Noy. 4 3 20.95 20 35.8 4.473 5.438 20.1 
22 22.16 11 38.7 11 3 16.15 20 19.7 
22 20.32 so.5 4438 4.762 1904 41 B30 HN Ole a ee Sue a0 
22 19.85 11 $7.9 Dec. 1 3 06.76 19 43.7 
22 20.74 1155.1 4.444 4.776 19.5 jee 02.87 9 19 27.0" 4,608 5.479 20.2 
22 22.94 11 44.4 21 2 59.90 19 13.0 
22 26.36 1126.4 “4.777 4.792" 19.7 31 2 58.05 1902.9 4.836 5.501 20.3 
22 30.86 11 01.8 Jan. 10 2 57.42 18 57.5 
22 36.30 1031.1 5.104 4.810 19.9 2078 2.58.05 {S570 20) 258189. WeSBI3H WO 
22 42.57 9 55.1 30 2 59.90 19 02.0 
22 49.51 914.4 5.398 4.831 20.0 | Feb. 9 3 02.92 194d 50 SUATG) 95,8260 2056 
22 57.01 8 29.6 19 3 07.01 19 25.2 
2304.95 — 741.4 5.638 4.855 20.1 | war 14 3 12.06 1942.6 5.813 5.569 20.8 
11 3 17.96 20 02.9 
63-64 21 3 24.61 2025.4 6.119 5.592 20.9 
: ; 025 31 3 31.89 20 49.4 
"y at eee a Aes pages ool ae Apreton 339.71. +21 14.3 6.373 5.616 21.0 
31. 0. 43.54 Spee 25.527) 51053 - 20.2 
ne 10 0 49.85 405.8 1965-66 
20 055.48 BGO) 51276) 5.076. 20.2 | Aug. 8 522.39 +2452.9 6.314 . 5.769 21.1 
30 100.34 5 22.9 18 5 28.99 25 00.3 
ly 10 104.32 553.4 4.995 5.094 20.1 98 +5 34.81 2506.4 6.064 5.796 21.0 
20 1 07.33 6 17.9 Sept. 7 5 39.75 5 11.6 
30 1 09.27 635.9 4.712 5.109° 20.0 17 5 43.70 9516.0. 5.783 \ 5.824 21.0 
ig. 9 1 10.09 6 47.1 27 5 46.55 25 19.8 
19 109.75 Grete 4568, 5,122. 19.9 | Oct. 7 5 48.21 9523.2 5.497 “5.881 20.9 
29 1 08.28 6 49.1 17 5 48.62 25 26.2 
pt. 8 105.75 640.4 4.262 5.134 19.8 27 5 47.76 9528.7 5.241 5.879 20.8 
18 1 02.34 6 26.3 Nov. 6 545.68 25 30.4 
28 058.28 608.1 4.160 5.145 19.7 16 542.48 D5 ait 252055) 5.907% 2027 
Ey) le 53.88 5 47.4 2 «5. 38.36 25 30.5 
18 049.48 Regge 47? 5.156. «19.7 | Dec 6  5'33.61 2528.4 4.965 5.936 20.7 
28 0 45.42 5 06.8 160 528785 25 24.6 
i 7» 42-02 450.9 4.300 5.166 19.8 26 = 5:23.56 9519.5 5.000 5.965 20.8 
ty = OBO 4 40.2 Jan. 5 5 18.99 25 13.4 
97 3RA2 ae ke 64596 - 5.477 ~~ 19:9 15 5 15.15 2507.0 5.156 5.994 20.8 
cc. 7 037.88 4 38.3 25 5 12.30 25 00.8 
17 0. 38.84 448.0 4.819 5.188 20.1 | Feb. 4 510.59 Bees A 53413) 46,023 224.0 
27, ~~: 40.97 5 04.8 14 510.10 4 51.2 
n 6 044.21 528.1 5.144 5.200 20.2 24 5 10.85 PHAR AQ OSS | 16.053. (21k 
16 048.47 587.5 Mar. 6 512.80 24 47.0 
Semen 6 371 «865.465 5.211 - 20.4 16 515.87 +2446.8 6.084 6.083 21.3 
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Taste III. Osculating values. 


Epoch M a) 2 1 q we oF f 
1960 Mar. 7 78°0 356°4 154°9 4°0 3.394 0.147 7293 i 
Sept. 23 100.9 358.6 154.4 4.0 3.404 0.148 7.98 ie 

1961 Apr. 11 121.8 2.9 153.0 4.1 3.416 0.150 8.06 OF 
Oct. 28 139.6 10.0 150.8 4.4 3.412 OStSom 8.12 0.. 

1962 May 16 1515 19.9 149.4 5.0 3.328 0.174. 8.09 0.. 
Dec, - 2 143.1 40.7 149.4 5.4 3.063 0.232 7.97 0.. 
1963 Jan. 11 13225 SoU 146.0 4.6 3.034 0.249 8.13 0. 
Feb. 20 111.8 95.9° 121.3 2.6 3.182 0.255 8.83 0. 

Apr. ot 68.4 233.8 25.6 3.9 3.875 0.244 11.6 0. 

May 11 22.9 302.2 12.6 7.0 4.672 0.320 18.0 0. 

June 20 9.8 327.4 12.6 7.0 4.958 0.385 22.9 0. 

July 30 6.4 338.0 1Sn0 6.2 5.053 0.393 24.0 0. 

Sept. 8 S155) 344.0 14.2 5.4 5.096 - 0.382 23%h Oy 

Oct. 18 5.6 348.3 14.3 4.7 5.12 . 0.367 23.0 0. 

Nov. 27 6.1 351.9 13.8 4.1 5.139 0.352 22.4 0. 

1964 Jan. 6 To S5Dee 12.9 3.6 5.154 0.339 21.8 0.. 
July 24 14.2 1355 125 223 Se205 0.292 20.0 0. 

1965 Feb. 9 23.0 33.5 345.6 1.9 5.273 0.268 19.3 0.. 
Aug. 28 31.9 46.4 335.9 1.9 5.322 0.256 19.2 OF 

1966 Mar. 16 41.0 52.1 332.4 1.9 5.383 0.250 19.2 1.1 
Oct) 2 50.4 54.5 $31.5 1.9 5.421 0.247 19.3 1.. 


be made afterwards, perhaps leading to an improved 
value for the mass of Jupiter. 
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